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Chapter 1 



Introduction 



1.1 Cosmology and gravitational waves 

It can be safely said that the advent of general relativity signaled the com- 
mencement of modern cosmology. The latter aims at a global description 
of the observable universe (i.e., the Universe). Based on observations at 
very large scales and some or other theory of gravitation cosmologists try 
to propose simple and plausible world models (i.e., models whose mathe- 
matical complexity is kept at a reasonable minimum and do not run into 
conflict with observational data). At present there are four main sources 
of observational data: the light from the faraway galaxies, whose redshifts 
seem to indicate that the Universe is expanding in an accelerated fashion; the 
abundance of light elements at cosmic scales and clusters thereof; the mass 
distribution of galaxies and clusters; and the anisotropies of the cosmic mi- 
crowave background radiation (CMB). These are usually interpreted in favor 
of the ACDM model, a flat Friedman- Robertson- Walker universe that began 
expanding from an initial phase of very high density and temperature (Big 
Bang), and today 70% of its energy would be in the form of a cosmological 
constant, A, and the rest in the form of cold (non-relativistic) matter. 

The future detection of gravitational waves (GWs henceforth) is expected 
to provide us with invaluable information about the instant of their decou- 
pling from other fields, i.e., about 1CT 43 seconds after the Big Bang. Likewise, 
GWs are a direct prediction of the Einstein's field equations whence it would 
entail a very strong evidence in favor of general relativity. Regrettably, the 
weak interaction of the gravitational field with matter -that renders GWs so 
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valuable from a cosmological viewpoint- hinders their detection. 

GWs can be found originated at single sources or in a background similar 
to the CMB. Of the first kind are, for instance, the GWs from supernovae 
burst and the periodic signals from spherically asymmetric neutron stars or 
quark stars; of the second kind are those GWs originated at the decay of 
cosmic strings and the relic GWs. 

This thesis is mainly concerned with relic GWs 1 , i.e., the GWs generated 
by parametric amplification of the quantum vacuum during the expansion 
of the Universe which form a background whose spectrum depends on the 
scale factor. The CMB is a remnant of the radiation that once dominated 
the expansion of the Universe and brings us information of the decoupling 
between radiation and non-relativistic matter (i.e., dust). In the same way 
the spectrum of the relic GWs can provide us with prime information about 
the evolution of the scale factor, i.e., about the history of the Universe. 

Even in the absence of any observed spectrum, some valuable information 
can be extracted from the existing bounds over it p. The presence of the 
relic GWs at the instant of the radiation-dust decoupling produces temper- 
ature fluctuations via the Sachs- Wolfe effect and, consequently, temperature 
anisotropies in the CMB. As these anisotriopies have been measured, a bound 
over the spectrum of relic GWs present at this instant follows. GWs passing 
by us transversely to our line of sight to a pulsar of well measured period 
will cause the arrival times of the pulses to shift. Many years of observation 
of the pulses arriving from a number of stable millisecond pulsars lead to 
another bound over the spectrum of the relic GWs. The theory of primordial 
nucleosyntesis predicts rather successfully the cosmic abundance of light el- 
ements. If at nucleosyntesis time, the contribution of relic GWs to the total 
energy were too large, the neutrons would have been more available resulting 
in an overproduction of helium. These restrictions over the GWs spectrum 
translates into constraints on the free parameters of the universe models (see 
top panel of Fig. II. 1|) . 

We assume throughout that GWs are negligibly damped by the cosmic 
medium; this is in keeping with common lore and well supported by different 
studies 



1 From chapter on, we refer to relic GWs just as GWs for the sake of simplicity. 
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1.2 Gravitational waves detection 

Because of insufficient technological means no direct detection of GWs has 
been possible so far. Nevertheless, some indirect evidences of their existence 
have been found. Binary systems (i.e., two stars orbiting their common center 
of mass) have an energy loss due to the emission of gravitational radiation 
that can be modelled by means of the general relativity. The observational 
data of the variation of the rotational period of the binary pulsar PS -R1913 + 
16 confirmed the predictions of the theory Recently, a similar but more 

precise indirect proof of the existence of GWs has been obtained from the 
data of the variation of the orbital period of the double pulsar J0737 — 3039 



The gravitational waves detection is a very topical subject and great effort 
is under way in that direction. There are two common types of detectors: 
resonant mass detectors and laser interferometers jH]- 

The first kind of detectors operates as the GWs produces a stress over 
a resonant bar (stretching or compressing it) and the vibrational amplitude 
or phase of the antenna experiments changes that can be measured. By 
making the antenna's quadrupole modes resonant at the wave's frequency, 
the detector keeps a "memory" of the excitation, allowing extra time to 
detect the signal. The sensitivity of this kind of antenna can be improved 
by increasing its mass (thus, increasing the force the GWs exert on it), by 
making it equally sensitive in all directions and polarizations and by lowering 
the thermal motion of molecules in the detector (thermal noise). Today there 
are seven cryogenic resonant-bar detectors in operation and three proposed 
resonant spheres. Analyzing the experimental results of the resonant bars 
EXPLORER and NAUTILUS, Astone et al. showed that the number of 
coincident detections is greatest when both of them are pointing into the 
center of our galaxy This conclusion, if not a direct detection of the 
GWs, can be considered a further indirect evidence of their existence. 

The laser interferometer detectors use test masses which are widely sep- 
arated and freely suspended (as pendulums reduce the effects of thermal 
noise); laser interferometry provides a means of sensing the motion of these 
masses produced as they interact with a gravitational wave. This technique 
is based on the Michelson interferometer and is particularly suited to the 
detection of gravitational waves as they have a quadrupole nature. Waves 
propagating perpendicular to the plane of the interferometer will result in 
one arm of the interferometer being increased in length while the other arm 
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is decreased and vice versa. The induced change in the length of the interfer- 
ometer arms results in a small change in the intensity of the light observed at 
the interferometer output 0. Nowadays there are four laser interferometers 
taking data at their level of sensitivity: TAMA300 jSj, GEO600 0, VIRGO 
[TO] and LIGO [llj (whose arms have a length of 300, 600, 3000 and 4000 
meters, respectively). 

Ground based detectors have a minimum frequency of detection bounded 
by seismic noise and atmospheric effects. Ground-based interferometers are 
also obviously limited in length to a few kilometers, restricting their coverage 
to events such as supernova core collapses and binary neutron star mergers. 
As a partial solution to these limitations, LISA - Laser Interferometer Space 
Antenna aimed at a launch in 2014 has been proposed by an Ameri- 
can/European team; it consists of an array of three drag free spacecraft at 
the vertices of an equilateral triangle of length of side 5 x 10 6 km. Proof 
masses inside the spacecraft (two in each spacecraft) form the end points 
of three separate but not independent interferometers. Each single two-arm 
Michelson type interferometer is formed from a vertex (actually consisting 
of the proof masses in a 'central' spacecraft), and the masses in two remote 
spacecraft. In the low-frequency band of LISA, sources are well known and 
signals are stable over long periods (many months to thousands of years). 
The primary sources for LISA are expected to be compact galactic binaries, 
supermassive black holes binaries, and extreme mass ratio inspirals. Other 
space based laser interferometers with a better sensitivity than LISA are be- 
ing proposed: the Advanced Laser Interferometer Antenna (ALIA) and the 
Big Bang Observer (BBO) [T3]. 

It is quite common, when tackling the subject of the G Ws detection, to 
define the dimensionless parameter Cl GW (f) as 



a Gw(f) = f 'l£i!V, (1 .d 

Pc 



where / denotes the frequency of the wave, p g is the energy density of the 
GWs and p c = 3H 2 / (8ir G). This parameter gives a measure of the intensity 
of the GWs signal and can be compared with the sensitivity of the GWs 
detectors (see bottom panel of Fig. II. lj) . 



1.2. GRAVITATIONAL WAVES DETECTION 



5 



±. -8 

51" -10 

o 

" - -12 
14 

-16 



_l 1 1 1 1 1 1 1 1 1 1 1 1 1 


i i i | i i | i i | i i | i i | i i | i i | i i | i i_ 
i 

/msec pulsar bound 


— 

: \cOBE bound 

" 


BBN bound- 

" 



-18-16-14 -12-10 -8 -6 -4 -2 2 4 6 S 10 
Log [f (Hz)] 




Figure 1.1: Top panel: Bounds over Clow from the CMB anisotropy (COBE), 
from the millisecond pulsar timing and from the Big Bang nucleosyntesis 
(BBN). Bottom panel: sensitivity versus frequency of LISA and ground-based 
detectors of first, second and third generations. The parameter ho is defined 
from Hq = ho x lOOkm/sec/Mpc and represents the existing experimental 
uncertainty in evaluating Hq. (Adapted from Figs. 1 and 2 of Ref. [Ti] ) 
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1.3 Outline of the thesis 

In chapter El we review the mechanism of amplification of the relic GWs and 
their power spectrum. In chapter EJ we obtain the relic GWs spectrum in a 
universe scenario consisting in a De Sitter stage of expansion followed by a 
radiation dominated stage and a dust dominated stage (three-stage model). 
Later, we obtain the spectrum in a four-stage model in which and additional 
era (dominated by a mixture of radiation and mini black holes) between the 
De Sitter and the radiation one. We compare both models and found bounds 
on the free parameters of the four-stage scenario. In chapter the spectrum 
of GWs is obtained in a scenario with an accelerated era of expansion domi- 
nated by dark energy right after the De Sitter-radiation-dust dominated eras. 
Although the spectrum is formally equal to those of the three-stage model of 
chapter El it evolves differently. The possibility of the existence of a second 
dust era if the accelerated era comes to an end in some future time is also 
considered. In chapter working in the De Sitter-radiation-dust-dark en- 
ergy scenario and assuming the GWs entropy is proportional to the number 
of GWs within the event horizon (being A the proportionality constant), we 
test the cosmological generalized second law of thermodynamics. An upper 
bound to the constant A is found. In chapter |U1 we leave the consideration of 
GWs to extend the study of the generalized second law of thermodynamics 
of the previous chapter to a universe dominated by dark energy. It turns 
out that the generalized second law is fulfilled in both phantom and non- 
phantom dark energy models provided the dark energy fluids have an entropy 
given by Gibbs' equation. Finally, the Summary discusses and sum up the 
main results of this work. Except where otherwise stated, we use the units 
such that H — G — c— ks — 1 • 

Before proceeding into the main body of the thesis it is sobering and 
expedient to recall the viewpoint of Willem de Sitter regarding cosmology 
research [To] : 

"It should not be forgotten that all this talk about the Universe 
involves a tremendous extrapolation, which is a very dangerous 
operation" . 



Chapter 2 

The amplification mechanism of 
gravitational waves 

2.1 The wave equation 

The existence of wave-like solutions of the linearized vacuum field equations, 
i.e., GWs, was first predicted by Einstein in 1916 when he realized of the 
propagation effects inherent in the gravitational field equations ^H]- The 
GWs equation to linear order was obtained by Lifshitz who considered per- 
turbations that not interfered with their own propagation (as they carry little 
energy and momentum) [Tj\. He proceeded by introducing a small pertur- 
bation in the background FRW metric, gij, 

9ij = 9ij + Kj(xi) = 0, 1, 2, 3), 

where \hij\ <C \gij\, and by evaluating the Einstein equations for the per- 
turbation to linear order in in a general coordinates system for which 

h 00 = h 0a = (a = 1,2,3). (2.1) 

Three kinds of solutions to the equations for the perturbation follow: scalar 
solutions, which represent density perturbations; vectorial solutions, which 
represent rotational perturbations; and tensorial solutions, which represent 
sourceless weak GWs. 

We focus our attention on the latter kind by considering a flat FRW 
universe with background metric 

ds 2 = -dt 2 + a(t) 2 [dr 2 + r 2 dn 2 ] = a^fl-drj 2 + dr 2 + r 2 dtl 2 ], (2.2) 
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where t and 77, the cosmic and conformal time, respectively, are related 
through the scale factor by a{rj)dr) = dt. 

Additionally to the conditions set by equation ()2.1|) . we impose the gauge 
conditions 

h lj .j = 0, (2.3) 

where the semi-colon denotes covariant derivative with respect to the back- 
ground metric ()2.2j) . From them, it is possible to show that the perturbation 
has two independent components only (i.e., the wave has just two polariza- 
tions) and that the wave equation simplifies to 

ha + + g^hj^s = 0, (2.4) 

a 

where the upper dot and the commas indicate partial derivatives respect to 
time and spatial coordinates, respectively. 
Introducing the ansatz 

hJ{t,x) = h{t)Gj(k,x), (2.5) 

where Gj 3 is a combination of plane- waves solutions exp(±2kx) that contain 
the two polarizations of the wave mentioned above and fulfill 

= -k 2 G?, G^ = G a a = 0, (2.6) 

into equation ()2.4|) leads to the evolution equation for the temporal part of 
the wave 

n ■ k 2 

h(k,t) + 3-h(k,t) + —h(k,t) = 0. (2.7) 
a a 2 

Or, in terms of the conformal time 

h"(k, 77) + 2-h'(k, 7]) + k 2 h(k, 7]) = 0, (2.8) 
a 

where the prime indicates derivative with respect to 77. 

In the equations of above we have used the comoving wave-number k = 
|k|, which is related with the wave-length and the frequency of the wave by 



2™ 

k = — — = aio. 



(2.9) 
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Equation ()2.8|) can be suitably simplified by using the ansatz h(k, 77) = 
fi(k,rj)/a(rj). This yields the so called Lifshitz's equation 

^" {v) + { k2 - a ~M) M = °- (2 - io) 

It parallels the time independent Schrodinger equation where the terms k 2 
and a" /a play the role of the energy and the potential, respectively. It can 
be also interpreted as the equation of an harmonic oscillator parametrically 
excited by the term a" /a. When k 2 3> i.e., for high frequency waves, ex- 
pression (j2.10j) becomes the equation of an harmonic oscillator whose solution 
is a free wave, and consequently the amplitude of h a/ 3(rj, x) decreases adia- 
batically as a -1 in an expanding universe. In the opposite regime, k 2 <C — , 
the general solution to (J2.1(Jj) is a linear combination of the pair of solutions 

/ii = a(r/), ^3 = 0(77) J drj a~ 2 (rj). (2.11) 

If the universe is expanding, [x\ will grow faster than /x 2 and will soon dom- 
inate. In this case, the amplitude of h a p(r],x) will remain constant as long 
as k 2 <C — . If at some future time this condition is no longer satisfied, the 
wave will have an amplitude larger than it would in accordance with the adia- 
batic behavior. This phenomenon is called "superadiabatic" or "parametric" 
amplification of gravitational waves |18j . 

GWs that fulfill the condition k 2 <C a" /a {k 2 ^> a" /a) are considered to 
be well inside (outside) of the Hubble radius [THj, i.e., their wave lengths are 
shorter (longer) than Xh = H" 1 ^), where H = a' /a 2 is the Hubble function, 
and 

k < 2naH(ri) (k > 2ixaH(r])) . (2.12) 

This can be traced to the fact that the ratio between (27raH(r])) 2 and the 
term a" /a is always constant for a perfect fluid with an equation of state of 
barotropic type, p — (7 — l)p, and tends to zero as soon as 7 — > 4/3. 

2.2 Gravitational waves creation 

In the previous subsection we have seen how the GWs may undergo "para- 
metric" amplification. This approach is based on the idea that the amplitude 
of each GW is enlarged during the expansion of the universe. Usually the 
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initial amplitude is assumed to be the vacuum one, A(k) ~ k [201 - ^ ne 
rationale behind this is the following. One assimilates the quantum zero- 
point fluctuations of vacuum with classical waves of certain amplitudes and 
arbitrary phases; consequently it is permissible to equalize u/2 with the en- 
ergy density of the gravitational waves, h 2 /X 2 times A 3 . Thereby the initial 
vacuum amplitude of GWs is given by h oc A(k) oc k [TS] . 

A different approach can be developed by realizing that Lifshitz's equation 
(j2.10|) is not invariant under conformal transformations of the metric 122] 

9a 9ij = tt 2 (t, x)g tj , (2.13) 

where Q 2 (t,x) is a real, continuous, finite and non vanishing function, except 
when a" = 0. With this in mind, the "parametric" amplification may be 
interpreted as spontaneous GWs creation due to the action of the expansion 
of the Universe over the initial vacuum. 

2.2.1 A simple example: the quantum pendulum 

The ideal pendulum, a weight hanging from a string of length I, pedagogically 
illustrates the phenomena of particle creation from the initial vacuum. The 
pendulum oscillates around its equilibrium position with frequency \fg~Jl- If 
the string is wound around a reel, the length of the string can vary as well 
as the oscillation frequency. 

Being the pendulum initially at the minimum energy configuration, we 
increase the length of the string by lowering the bob, in a time interval 
t = tf — ti, from an initial value Zj to a final value If. Classically, the 
pendulum (initially at rest) will be also at rest in the final state and the 
change of potential energy will equal the friction work. 

From a quantum mechanic point of view, the initial minimum energy 
configuration of the system corresponds to the state with energy Ei = Ui/2 
and number of quanta iVj = 0. If the time spent in lowering the bob is much 
longer than the initial and final periods of oscillation, i.e., r ^> Tf = 2n / 'ujf 
(the adiabatic case), the bob undergoes more than one complete oscillation 
during the process. In this case the final energy is Ef — Uf/2 and the number 
of quanta in the final state is Nf = 0. The potential energy obtained is partly 
dissipated by the reel, as in the classical approach, but also by the change of 
the energy from the initial to the final state. 

By contrast, if the time spent in lowering the bob is much shorter than 
the initial and final periods of oscillation, i.e., r <C Ti = 2ir/u)i, things fare 
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differently. The final energy of the oscillator will be half the energy of the 
initial state Ef = Ei/2 and the final vacuum energy will be uJf/2 (see Ref [I] 
for details). Consequently, the final number of quanta Nf = wm differ 
from zero. Note that these quanta are created because of the sudden change 
in the vacuum state. 

In the case of the GWs, we must replace the pendulum equation for 
equation (|2.10j) . The role of the variable length I is now played by the term 
a" (77) / a(rj) . If a"{rj)/a(rj) varies slowly compared with the wave number k, 
the initial vacuum state will evolve into the final vacuum state and no GW 
will be generated. But if a"(rj)/a(ij) evolves fast enough, the final state will 
not be the vacuum one and GWs will be spontaneously generated. 

2.2.2 Bogoliubov coefficients 

In Minkowski spacetime, a scalar field <fi of mass m obeys the Klein-Gordon 
equation 

(rj ii d i d j +m 2 )<f> = 0, (2.14) 
where rfi represents the metric. The field <fi can be expanded as 

(f>(x) = J2l A ^( x ) + 4X(z)], (2-15) 

k 

in terms of the family of modes 

Z/ k = l_ e i(k-x-U) ; ( 216 ) 

which are positive-frequency defined with respect to t 

<9z/ k 



dt 



-ikv k (k > 0). (2.17) 



This set of modes forms a privileged family: it defines creation and annihi- 
lation operators A k and A^, respectively, representing real particles for all 
inertial observers, as the quantum vacuum defined by the modes is invariant 
under Poincare transformations. 

In curved spacetime, the equation obeyed by the scalar field reads 

gv^i.j + (m 2 + £R) (f) = (2.18) 
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where £ is a dimensionless parameter and R is the Ricci scalar. The term £i? 
accounts for the coupling between the scalar and the gravitational field. 

When m = and £ = 1/6, equation (J2.18)) is conformally invariant. 
In this case if the spacetime considered can be transformed conformally in 
the Minkowski one, it is possible to define a new field that obeys equation 
(|2.14j) and privileged modes can be chosen in order to define a vacuum state. 
However, generally, there is no privileged family of modes as the Poincare 
group is no longer a symmetry. 

In general, given a complete collection of modes v m (x) the field can be 
expanded as 

:) = Y}4n* m (x)+Altf n (x)], (2.19) 



where A m and A^ are, respectively, the annihilation and creation operators 
associated to the quantum vacuum state of the family, |0). The field can also 
be expanded as 



[x 



J2[^nV n (x) + Alv* n (x)], (2.20) 



in terms of a different complete collection of modes, V n (x), with a different 
vacuum state, |0). 

The Bogoliubov transformation 



u n ix) 



^ \oL nm Vm{x) + f3 nm v* m {x)} , (2.21) 



where a nm and (3 nm are the coefficients of Bogoliubov, relates both families 
of modes. From equations ()2.19|) - (j2.21|) . it is readily seen that 

An = ^(a nm Z n + /3 nm Zl), (2.22) 

n 

A-n = ^ ^ (&nmA m + (3nm^\n) ■ (2.23) 



As the creation/ annihilation operators satisfy the commutation relation 
[y4 m ,y4j„] = [y4 n ,y4^] = 1, the Bogoliubov coefficients must fulfill 

a nm \ 2 - \P nm \ 2 = 1. (2.24) 
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The operator number of particles of the first family, N m = A^ m A m , acts 
upon the quantum vacuum state of the second family according to 

(0|Ayo> = (0\AlA m \0) = \Pnm\ 2 . (2.25) 

n 

Therefore, the quantum vacuum state of the second family contains J2 n I Aim 1 2 
particles of the first family. An entirely similar relation holds for the first 
family vacuum state and the second family number of particles operator. 

The GWs equation ()2.8|) formally coincides with the equation of a mass- 
less decoupled field in a spatially flat FRW spacetime. Its solutions are 

fcyfax) = J i^A {k) hf{k ) x)+A\ k) hf j ) \k ) x)) d 3 k. (2.26) 
The modes hff (k, x) can be written as 

where ey(k) contains both polarizations of the wave, and the functions P(k){v) 
are solutions to (|2.1(J|) . As the family of modes that quantify the field is 
complete, the modes h\j\k,x) must be orthogonal, i.e., 

h%\k,x)hf; V '(k,x) - hf(k,x)h^\k,x)) d n - l x = 5(k - k'), 



(2.28) 

and, therefore, p^) obeys the further condition 

M(fe)M(fe) - ^(k)^(fc) = L ( 2 - 29 ) 

If the spatially flat FRW universe contains a perfect fluid with equation 
of state 

V = (7 - l)p, (2-30) 
with 7 constant, the Einstein equations 



a'\ 2 Stt 
~7i 



p, (2.31) 



p' + 3-(p + p) = 0, (2.32) 
a 
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djHi \ ( I x 



can be solved to 



where a iy Hi and 77^ are, respectively, the initial values of the scale factor, the 
Hubble function and the conformal time and / = 2/(37 — 2). For simplicity, 
we define r)i — r\ — % + //(a^if;) and Lifshitz's equation reduces to the Bessel 
equation 



k 



2 Ki-K 



Kv) = 0, (2.34) 



(Vi) 2 

whose solutions can be written in terms of Hankel functions as 

M(V) = ^-^x] 12 \D U H^\ ( Xl ) + D 2l H {2 \ ( Xl )] , (2.35) 
2 L 1 2 1 2 

where x\ = krji, and Du and D 2 i are integration constants. 

Henceforward, we shall work in the Heisenberg picture, the quantum 
states are functions of time meanwhile the operators are time independent. 
Initially the universe will be in some state from which the number of GWs 
operator will be defined. As the scale factor evolves, the quantum state 
evolves too and the constant operator acts upon this state which is no longer 
the initial one. 

In this scenario, it seems problematic to choose the appropriate vacuum 
state or to define real particles. A possible solution to this problem rests on 
the concept of adiabatic vacuum. 

2.2.3 The adiabatic vacuum approximation 

This approximation rests on the idea that the creation of particles by a slow 
change in the initial state is minimal [22j. 
Solutions to the equation 

x " + uj 2 (r ] )x = 0, (2.36) 

with uo G C°° are expressed as a linear combination of basis functions which 
can be approximated by 

X ± (r ] ) = y=y ex P (=F* /" W(rf)drf) + 0(oj~ n ), (2.37) 
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with 

W 2 ( V ) = [I + X 2(v)u- 2 + X*(v)u~ 4 + •••] , (2.38) 

where the Xn(v) are functions of uj{rj) and its derivatives up to the derivative 
of order n, which is bounded as u — ► oo [23 . If we introduce the adiabatic 
parameter by replacing Qu(r]) for u(rj) (at the end we can let 6 = 1), the 
adiabatic behavior (slow expansion limit) of equation ()2.36|) can be examined 
when — > oo, and, in this limit, the above expression can be considered a 
power-series expansion in _1 to order N. 

Equations (j2.1()j) and ()2.36j) share the same form. Thus, we can define 
the adiabatic vacuum of order iV by taking [i in equation ()2.27|) as an exact 
solution to ()2.10|) whose WKB approximation is precisely x + m f)2.37|) . To 
this order, the operators A(k) in ()2.26j) correspond exactly to physical particles 
when — > oo or when k — > oo. 

(k) 

Note that in this limit, the modes h\j tend to the positive-frequency 
Minkowski modes. Hence, the constants in equation ()2.35|) turn out to be 
Du = and D 2 i = 1, thereby 

fit = {V^me^k- I l 2 x) /2 H\ 2 \ ( Xl ). (2.39) 

1 2 

2.2.4 The creation mechanism 

Let us assume that the spatially flat FRW universe considered above goes 
through a succession of stages, and that, in a generic stage r, the Universe is 
dominated by a perfect fluid of barotropic index 7 r . The transitions between 
the different stages can be considered instantaneous or, more accurately, 
much shorter than the duration of each stage. 

The scale factor, in this scenario, is given by equation ()2.33|) with a dif- 
ferent l r in each era and the subindex i denoting the beginning of the r-stage 
(i.e., the end of the (r — l)-stage). Note that the scale factor must be contin- 
uous at each transition rji as no discontinuities must be present in the FRW 
metric. In each stage the solution to Lifshitz's equation ()2.10j) is (|2.39j) with 
/ replaced by l r . 

The modes of two consecutive stages (namely the (r — l)-stage and the 
r-stage) are related by a Bogoliubov transformation 

Mr-M = atnUvi) + MM)- ( 2 - 4 o) 

Since the functions Hi r {ji) an d (v) must be also continuous at rji (fur- 
ther details about the continuity of /i and the validity of the sudden transition 
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approximation are given in the Appendix), it is possible to obtain the ex- 
pression of the Bogoliubov coefficients in terms of known quantities |2*4*1 12*5] 

a i = ilPir^ivdvUVi) ~ Vir-Ml4 r '(Vi)}, (2-41) 

The adiabatic approximation sets a bound over the GWs that can be 
created in this transition. The limit of slow expansion, u —>■ oo, can be 
relaxed by assuming that there is no creation of waves for u(rji) ^> w^. This 
bound can be fixed as the frequency associated to the characteristic time in 
which takes place the transition. Modes whose period at the instant rji are 
shorter than this transition time experiment an adiabatic expansion thereby 
they still represent the quantum vacuum state, i.e., for them, a, = 1 and 
(3i = 0. Meanwhile modes whose period at the transition rji is larger than this 
characteristic time do not represent the quantum vacuum state any longer, 
and f3i ^ 0. This characteristic time is usually chosen as the inverse of the 
Hubble parameter at the transition, Hj~ . Then, the adiabatic bound for the 
frequencies at rj > rji reads 

Ui(ri) = 2ir-^H h (2.43) 
a(rj) 

where we have taken into account the redshift. 

Assuming that each created GW has an energy 2u(rj) (the factor 2 
comes from the two polarization states), it is possible to express the energy 
density of GWs created at the transition rji with frequencies in the range 
[uj(r}),u(r]) + duj{rj)} as 

|A| 2 = P(uj( V ))du( V ), (2.44) 

where P(u(t))) = (cu 3 (rj) / V 2 ) |/3,| 2 denotes the power spectrum. As the energy 
density is a locally defined quantity, p g loses its meaning for metric pertur- 
bations whose wave length A = 2n /uj{rj) exceed the Hubble radius H~ l {rj). 
Integrating over the frequency, we get the total energy density 



dpgiv) = 2a; (77) 



u (V) 
2tt 2 



duj(rj) 



Pg(v) = / P(u(rj))du){ri). 

J2nH(ri) 



(2.45) 
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The adiabatic cutoff amounts to a renormalization of the quantum field the- 
ory. This upper bound ensures that the energy density does not diverge at 
high frequencies. The adiabatic bound plays the same role that the extrac- 
tion of the energy of the vacuum in the quantum field theories in Minkowski 
spacetime. 

To relate the \i\ functions of two non-consecutive stages, we must make 
use of the total Bogoliubov coefficients which can be found recursively from 

a T , = OLidT^ + PiPr^ v (2-46) 



@Ti — A a T {i _i) + a i /^Tfi-i) , (2.47) 

where the sub index (z — 1) denotes the previous transition to the i-th one 

Equations (|2.4fij) and ()2.47j) readily follow from the Bogoliubov transfor- 
mation 

(&)-(« «)(*)■ (248) 

Thus, the modes of the r — n stage, /i; r _ n , are related with 



( M r _ n \ ( Qfi-(n-l) A-(n-l) \ ( Qti_i A-l \ ( a i & \ ( Mr \ 

(2.49) 

(note there is n — 1 transitions between both eras). After simplifying we get 



/i; r _„ \ = ( a Tl fa \ ( Mr 
/£_„ J V fti a h J V M* r 
We will employ Eqs. (|2.4fij) and (|2.47jl in the following chapters. 



(2.50) 
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Chapter 3 

Quantum mini black holes and 
the gravitational waves 
spectrum 



In this chapter we consider in detail the power spectrum of relic GWs assum- 
ing the Universe went through an expansion era dominated by radiation and 
mini black holes (MBHs) intermediate between the conventional inflationary 
(De Sitter) era and the radiation dominated era |26| . The existence of that 
era is justified in subsection 13.21 In subsection 13 . 1 1 we recall the derivation of 
the power spectrum of the conventional three-stage scenario. 

3.1 Three-stage spatially flat FRW scenario 

Here we apply the fundamental equations found in the previous chapter to a 
spatially flat FRW scenario consisting in an initial De Sitter stage, followed 
by a radiation dominated era and a dust era that includes the present time. 
The power spectrum predicted for this scenario is well known in the literature, 
see Refs. [2I| and [211-123 ■ The scale factor is given by 



The initial state is the vacuum associated with the modes of the inflation- 
ary stage fii(fj), which are a solution to Lifshitz's equation ()2.1()j) compatible 




(-00 < rj < 771 < 0), 
1i ) (vi <V < V2), 
^- (V2 < V < Vo)- 



(3.1) 
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with the condition ()2.29j) . Taking into account that I = — 1 in this era (De 
Sitter), the modes read 

m = (^/2)e i ^k- 1 l 2 x l / 2 H% /2 (x) 1 (3.2) 

(2) 

where x = kr], ipj is an arbitrary constant phase and H_L 2 {x) is the Hankel 
function of order —3/2. The proper modes of the radiation era (I = 1) are 

Ai fl = {^/2y^k- x ' 2 xfHf/ 2 {x R l (3.3) 

where xr = k(rj — 2rji) and i/jr is again a constant phase. 

The "parametric" amplification in this first transition was first developed 
by Grishchuk [21] and Starobinsky j2H]- Following the quantum approach, 
the two families of modes are related by 

fjn(ri) = ai^ R {ri) + ft/i^fa). (3.4) 

From equations (j2.41j) and (J2.42|) . we obtain 

i 1 1 

where we have neglected an irrelevant phase. Considering the adiabatic quan- 
tum approximation, the modes whose frequency at the transition are larger 
than the characteristic time scale (H^ 1 ) get exponentially suppressed. Thus, 
the coefficients are a± — 1 and /?i = for GWs with k > 2iraiHi and (J3.5|) 
when k < 27ra 1 H 1 . 

In the dust era (rj > r/ 2 and I = 2) the solution for the modes is 

fi D = (V^me^k-^x^H^xo), (3.6) 

where xd = k (r] + i]2 — 4r^i) and it is related to the radiation ones by 

Vr(v) = ®2Vd(v) + foUoiv)- (3-7) 

Similarly one obtains 

0(2 = - i ( 1 + 2k( V2 - 2 Vl ) ~ 8(k(r ]2 -2r ]l )) 2 ) ' ^ = 8 (k^ - 2 Vl )) 2 ' 

(3.8) 
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for k < 27ia(r]2)H2 and 0L2 = 1, /?2 = for k > 2-nairi2)H2 where H2 is the 
Hubble function evaluated at the radiation-matter transition 772. 

In order to relate the modes of the inflationary era to the modes of the 
dust era, we make use of the total Bogoliubov coefficients ar r 2 and /?rr2, given 
by equations (|2.46|) and (|2.47j) . For k > 2^a\Hi, we find that 0^2 = 1 and 
@Tr2 — 0; in the range 2 r na\H\ > k > 2HCL2H2, the coefficients are aj r2 = «i 
and (3tt2 — Pi, and finally for k < 2na(ri2)H2 we obtain that 

2(772 -2rn) +771 
" "8^(772 -2r /l ) 2 ' ( } 

Thus the number of GWs at the present time, 770, created from the initial 
vacuum state is (Nj) = \PTr2f ~ ^ 6 (vo) f° r ^(Vo) < 2ir(a2/ao)H2, ^~ 4 ( ? 7o) 
for 2ir(ai/a )Hi > ^(770) > 2ti (02 / ao) H2 and zero for 0^(770) > 2n(ai/ao)Hi, 
where we have used the present value of the frequency, ^(770) = k/ao. 

The present power spectrum of GWs predicted by this model is 

( (uj(rj ) > 2ir(a 1 /a )H 1 ) , 

P(w) ~ < u- l (r] ) (27i(a 2 /ao)H2 < u( m ) < 2 7 r(a 1 /a )# 1 ) , (3.10) 
[ ^- 3 (t7o) (2ttH < u(r) ) < 2tt (a 2 / a ) H 2 ) ■ 

The spectrum (|3.10j) is shown in figure 13.11 tagged as I = 1 for different 
choices of Hi, ai/ao and 02/00- 

Below we compare the predictions of the three-stage model of above, for 
the frequency ranges and powers, with the four-stage model that includes 
the "MBHs+rad" era. 



3.2 GWs in a FRW universe with an era of 
mini black holes and radiation 

3.2.1 Mini black holes in the very early Universe 

As is well known, MBHs can be created by quantum tunnelling from the hot 
radiation with a rate of nucleation given by jHO] 

r(T) = 0.87T(^)^e- 1 / 16 ^ 2 , (3.11) 

where T is the temperature of the radiation, \i is a parameter close to the 
Planck mass, and 6 is a numerical factor which depends on the number of 
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spin fields accessible to the system (i.e., 6 = 3.083 for the standard model, 
3.656 for the supersymmetric standard model, 3.000 for the supersymmetric 
SU(5) and 0.283 for the SU(5)). 

Thus, the number of MBHs per unit comoving volume at time t is 



t 

3/ 



n{t)~— J^a\t')r{T{t'))dt', (3.12) 

where t* is the time of formation of a MBH which would have just evaporated 
at time t. And the density of MBHs is 

i r* 

Pbh = - a 3 (t')M(t,t f )T(T(t'))dt\ (3.13) 
a Jt* 

where M(t, t') is the mass a MBH would have by time t if formed at time t'. 

The MBHs nucleated by quantum tunnelling will be strongly peaked 
around an initial mass given by 



M. 



BH\t 



form 



8nT 



(3.14) 



t f 



form 



where tf orm is the instant of the MBH formation and T is the temperature 
of the radiation. Assuming that during the formation of the MBHs the 
expansion of the Universe is dominated by the hot radiation 

P = § , P = ^9*T\ (3.15) 

here g* is the number of relativistic particles species, and using the Friedman 

equations, one has that T oc t -1 / 2 . Consequently, the initial mass of MBHs 
i li 

will increase as tj orm . The evolution of the mass of the MBH from this point 
onward will depend on its interactions with the surrounding radiation and 
with the other MBHs [2]. 

If the total number of MBHs in a comoving volume is constant, the rela- 
tive velocity between two neighbor MBHs due to the expansion is 



~ 1/3 H, (3.16) 



where n is the MBHs number density. The velocity required for a MBH to 
escape the gravitational attraction of its neighbor is 



v esc = V Amn l / Z , 



(3.17) 
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where m = psu/n is the mass of typical MBH. Given that during this period 
H = a/ (87r/3)(p ra rf + Pbh), it is straightforward to conclude that v exp > v esc . 
Therefore, the collisions between neighbor MBHs are not frequent enough to 
significantly affect the mass spectrum of MBHs. 

Even though the MBHs form at the same temperature as their surround- 
ings, one might imagine that cooling of the radiation due to Hubble expan- 
sion would cause the MBH to begin to evaporate freely immediately after 
formation. Therefore, their mass would evolve according to 



{dM/dt) free = -W(3M 2 ), (3-18) 

where g** is the number of particle species for the black hole to evaporate 
into. 

But the MBH evaporation can be delayed if it begins absorbing radi- 
ation after its formation so that its mass increases. Comparing the aver- 
age time of interaction between the MBHs and the surrounding radiation, 
T int = l/(367rm/3B#), with the characteristic time of expansion of the Uni- 
verse (the inverse of the Hubble factor, t exp = H' 1 ), it is straightforward 
to demonstrate that, although these interactions are not frequent enough to 
consider an evolution different from (j3.18j) for all the particle models of inter- 
est, the energy density of the MBHs can become comparable or even exceed 
that of the radiation for sufficiently high temperatures |31j . 

Further, it is also natural to assume that the MBHs are surrounded by 
an atmosphere of particles in quasi thermal equilibrium with them The 
MBH emits quanta in a perfectly thermal manner and these quanta might 
create a thermal atmosphere surrounding tightly the MBH. The vast ma- 
jority of the quanta emitted are prevented from escape the MBH gravita- 
tional potential whether because they have a large angular momentum or 
because, having an adequately small angular momentum, they tend to have 
such large wave lengths that when trying to escape they scatter of the MBH 
spacetime curvature and are driven back towards the horizon. Thus the 
absorption and emission of particles from the atmosphere to the MBH and 
vice versa prevents the MBH to evaporate freely. Therefore, we have that 
\(dM/dt) atm \^\(dM/dt) free . 

Even assuming the MBHs begin to evaporate freely since they are cre- 
ated, except in the SU(5) model, pbh will be comparable to the radiation 
density in a time span of two to one hundred times the Planck time from 
the instant the nucleation starts jHH|33]. It is reasonable to expect that, at 
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this point, a steady state would be achieved where the total energy density 
is shared between the black holes and the radiation whence p = Pbh + Pr 
and, consequently, the total pressure is 

P = Pbh + Pr = (l - l)p, (3.19) 

where the constant 7 lies in the interval 1 < 7 < 4/3. If the density of MBHs 
is large enough to dominate the expansion of the Universe, then 7 ~ 1. In 
the opposite case, the Universe expansion is dominated by the radiation, 
7 ~ 4/3. From the Einstein equations and (j3.19j) one finds that 1 < I < 2 
during the "MBHs+rad" era. The MBHs eventually evaporate in relativistic 
particles after some time span which can be fairly large if MBHs have a 
thermal atmosphere. 

Interestingly, the evaporation of MBHs when the age of the Universe 
was about 100 Planck times may explain why the cosmic baryon-number to 
photon ratio is of the order of 10~ 9 . Lindley argued that that small figure 
can be obtained if the Universe's expansion was dominated by black holes of 
a few hundred Planck mass at the mentioned epoch [S3]. This looks feasible 
in the scenario contemplated here. 



3.2.2 Power spectrum of the four-stage scenario 

In this subsection we obtain the power spectrum of the GWs assuming the 
following eras in succession: an initial De Sitter era, an era dominated by 
MBHs and radiation, the conventional radiation-dominated era and the dust 



era. 



The scale factor of this four-stage scenario is 



1 

~H 1V 



a{<n) 



where tjbh 
n C+i) 



(—00 < r] < rji < 0), (De Sitter era) 
(771 < 77 < 972), ("MBHs+rad" era) 
(772 < 77 < 773), (radiation era) 
[Vd] 2 (m < V < Vo), (dust era) 

(3.20) 

77- (/ + 1)771, r} R = 77 + ^^772-^^771, t] D = 77 + 773 + 2^^7/2- 



i i Hl(-77i) i + 1 
Hi(-»?i) ,+1 lR 

2H 1 (- Vl y+ i VDa 



f^Vi, VR2 = [V2 - (I + l)77i] /I and 77^3 = 2 
the previous section, the sudden transition approximation is assumed 



, (i-O 

7?3 + ??2 



(*+l) 



As in 
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The shape of fi(r)) can be found by solving Lifshitz's equation (|2.1()jl in 
each era. For the De Sitter era, fi(r]) is given by (|3.2|) as above. For the 
"MBHs +rad " era the solution of (EHUJ) is 

llBH = (V^W^k-^X^H^ixBH), (3.21) 

where xbh = kr/BH, and it is related to the modes of inflation by the Bo- 
goliubov transformation ()3.4j) with fM R replaced by hbh- By evaluating the 
Bogoliubov coefficients from ()2.41j) and ()2.42j) . we obtain 



7T 1 



H^likrj^H^i-lkr,,) - H^-lk^H^k^) 



(3.22) 



# = ~Z3(fc77i) H^lik^H^i-lkv^-Hp^-lk^H^ikm) . (3.23) 



r(2), 



r(2) 



(2) 



r(2), 



In the small argument limit (i.e., x <C 1) the Hankel functions can be 
approximated by 




(2m- 1)! 



7TX ' (m — I)! 

Thus, the coefficients dominant term is 

l 2 2 l 



T f-(^2ix)- m+ 2e 



(3.24) 



[-l V 



(3.25) 



when /c < 27ra 1 if 1 and = 1, /3[ = when /c > 2-na\H\. 

The solution for the radiation era is again ()3.3|) with x# = /c?7j?. The 
coefficients that relate (|3.3|) with \i B h are 



«0 



1 JttIx R2 



(J- -i]H\ 2 \{lx R2 )-H ( fUlx m ) 

\X R2 J 1 2 '+2 



o«Zfl2 



(3.26) 



1 j7dXR2 



(— + ^ ff. (2 \ (/x R2 ) - ilx m ) 

\XR2 J l S '+2 



-*Zfl2 



(3.27) 



when < 27ra2-£/2 and ai = 1, (3\ = when > litaiH.^ 
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The modes of the dust era are given by (jH.fijl . with xd = kr)z>, and they 
are related to the modes of the radiation era by the coefficients 

i 1 \ . 1 



when /c < 2Tra^H^ and 03 = 1, = when k > 

The total coefficients relating the initial vacuum state with the modes of 
the radiation state can be evaluated from (|2.46|) and (|2.47|) . In this case, we 
get 

a l Tr2 = a l 2 a[ + P2PI = $Wi + aftPl (3.29) 
and, consequently, the total coefficients dominant term is 

1,0 (k> 2ira 1 H 1 ) , 

/?' ~ I <*{, (3[ {2na 1 H 1 >k> 2na 2 H 2 ) , ( 3 . 30 ) 



a: 



Tr2> M7V2 



^i) i+1 (^2) 



where we have made use of (|3.24j) in order to approximate the coefficients a l 2 
and (3 l 2 . 

Finally, the total coefficients relating the inflationary modes with the 
modes of the dust era evaluated from 



a 



Tr3 



a 3 a l Tr2 + (3lf3 l Tr2 , (3 l Tr3 = (3 s a l Tr2 + a\(5 l Tr2 , (3.31) 



are found to be 

a Tr2i /3 l Tr2 (k > 2ira 3 H : 

(2l 2 -3l+l) 



We are now in position to calculate the current spectrum of GWs. Taking 
into account that 

771 = -(oi-Hi) -1 , 



1/1 



VR2=[f 1 ) (ai/fi)- 1 , (3.33) 
VD3 = 2f 2 (f) V ' {a x H x r\ 
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and ui = k/ao, the GWs power spectrum can be written as 



a 



P(u) ~ < 



i a ""(2i 2 -3f+i) 2 faA 41 (a^ 2 - 2/l 



2tt^Hx >uj> 2tt^H 2 

an 1 a z 



64?r 2 



16^ 2 ' 



do J \ a2 

4Z+4 / \ 2 



H u -(41-3) ( 2tt^H 2 >uj> 2ir^fr 

\ ao z an 



an ' 



(3.34) 



Comparing the power of u in (|3.10jl and (|3.34|) for uj < 2it^Hi, we con- 
clude that the four-stage scenario leads to a higher number of GWs created 
at low frequencies than the three-stage scenario. This fact can be explained 
intuitively with the classical amplification approach. In the three-stage sce- 
nario the GWs are parametrically amplified as long as k 2 < a" {rj) / a{rj) . For 
7] = T]x, a"(j])/a(j]) vanishes and there is no further amplification. On the 
other hand, in the four-stage scenario the GWs with u < 2it—Hi are am- 
plified until the instant rji by identical term a"(rj)/a(rj) = 2/rf than in the 
three-stage scenario and from r]\ to r] 2 by the term 1(1 — 1)/ rf BH . Consequently 
they have a larger amplitude in the radiation era. 

Figure 13.11 shows the spectrum (|3.34|) for I = 2 and I — 1.1. As is 
apparent the four-stage scenario gives rise to a much lower power spectrum 
than the three-stage scenario assuming that in each case the spectrum has 
the maximum value allowed by the CMB bound. The higher the MBHs 
contribution to the energy density, the lower the final power spectrum. 

Thus, if LISA fails to detect a spectrum at the level expected by the 
three-stage model, rather than signaling that the recycle model of Khoury et 
al. jSZ] 1 should supersede the standard big-bang inflationary model it may 
indicate a "MBHs+radiation" era between inflation and radiation dominance 
truly took place. Likewise, once the spectrum is successfully measured we 

1 Khoury et al. have proposed an alternative brane- world cosmological scenario which 
addresses the cosmological horizon, flatness and monopole problems an generates a nearly 
scale-invariant spectrum of density perturbations without invoking an inflationary period. 
In that model, the spectrum of GWs is strongly blue in comparison with these of the 
standard big-bang inflationary model. 
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will be able to learn from it the proportion of MBHs and radiation in the 
mixture phase. 

In this four-stage cosmological scenario, the Hubble function H{rj) de- 
creases monotonically, while the energy density of the GWs for 77 > 773 can 
be approximated by p g {f]) ~ H~ Al+2 (r]) thereby it increases with expansion 
|21]. Obviously this scenario will break down before p g (rj) becomes com- 
parable to the energy density of matter and/or radiation since from that 
moment on the linear approximation on which our approach is based ceases 
to be valid. In the two next subsections constraints on the parameters of the 
model are imposed so that this does not happen. 



3.2.3 Free parameters of the four-stage scenario 

At this stage it is expedient to evaluate the parameters occurring in (|3.34j) . 
The redshift relating the present value of the scale factor with the scale 
factor at the transition radiation-dust, may be taken as 10 4 [HE]. The Hubble 
factor Hi is connected to the energy density at the inflationary era by 

* = l^Hl (3.35) 

where we have restored momentarily the fundamental constants. In any rea- 
sonable model the energy density at that time must be larger than the nuclear 
density (~ 10 35 erg /cm 3 ) and lower than the Planck density (~ 10 115 erg /cm 3 ) 
therefore 

lOV 1 <H X < lO 43 ^ 1 . (3.36) 

Using the expression for the scale factor at the "MBHs+rad" era in terms 
of the proper time, one obtains 

n / 111 \ '/(l+l) 

^= l + H^r) , (3.37) 



where r is the time span of the "MBHs+rad" era which depends on the 
evaporation history of the MBHs. 

However, the "MBHs+rad" era span r should be longer than the dura- 
tion of the transition at rji (as the transition is assumed instantaneous) in 
calculating the spectrum of GWs. To evaluate the adiabatic vacuum cutoff 
for the frequency we have considered that the transition between whatever 
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two successive stages has a duration of the same order as the inverse of the 
Hubble factor. This places the additional constraint 

r > Hi 1 . (3.38) 

Finally, ^ can be evaluated from the evolution of the Hubble factor until 



„ N V2 / \ {l-W /„ \ 2 



the present time 



The current value of the Hubble factor is estimated to be 2.24 x 10~ 18 s _1 
and 

n \ 2 (n \ 1/2 / 7-1-1 \ (i-O/C+i) tt 

.|)-(S + |; «-°) 

The only free parameters considered here are I, r and H±, with the restric- 
tions 1 < I < 2, (|3.36|) and ()3.38|) . The two first free parameters depend on 
the assumption made on the MBHs, although it is possible to obtain rigorous 
constraints on Hi and r from the density of the GWs. 



3.2.4 Restrictions on the "MBHs+rad" era from the 
cosmic microwave background 

It is obvious that p g cannot be arbitrarily large, in fact the GWs are seen as 
linear perturbations of the metric. The linear approximation holds only for 
p g {rj) p(rj), p{rj) being the total energy density of the Universe. Several 
observational data place constraints on p g . The regularity of the pulses of 
stable millisecond pulsars sets a constraint at frequencies of order 10~ 8 Hz 
|4Uj . Likewise, there is a certain maximum value for p g compatible with 
the primordial nucleosynthesis scenario. But the most severe constraints 
come from the high isotropy degree of the CMB. We will focus on the lat- 
ter constraint. The observed thermal fluctuations are usually analyzed by 
decomposing them into spherical harmonics 

ST 

^T = ^a, m F, m (M), (3-41) 

where ai m are expansion coefficients and 9 and (f) are spherical polar angles 
on the sky. Defining the power spectrum by Ce = (\ae m \ 2 ), it is conventional 
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log 13 (m.is) log 15 (">l£) 

Figure 3.1: GWs spectrum for an expanding universe with a "MBHs+rad" 
era for certain values of /, r and H±. The spectrum predicted for the three- 
stage model of the previous section is plotted for comparison, 1 = 1. Pa- 
rameters r and H 1 are chosen assuming each spectrum has the maximum 
value allowed by the CMB anisotropy data at the frequency uo = 2ttH = 
2.24 x 10 _18 s _1 . In the bottom-right panel the power spectrum with I — 2 is 
rule out as it yields a CMB anisotropy larger than observed. 
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Figure 3.2: The CMB power spectrum from the WMAP satellite jlj- The 
error bars on this plot are 1-a and the solid line represents the best-fit cos- 
mological model 13D: . 



to plot £(£+l)Ce versus £. Figure l3~2l shows the first published results of the 
WMAP experiment regarding the CMB temperature anisotropies as well as 
the ACDM model (solid line) that fits the data rather well. 

Metric perturbations with frequencies between 10 -16 and 10~ 18 Hz at the 
last scattering surface can produce thermal fluctuations in the CMB due to 
the Sachs- Wolfe effect (323- These fluctuations cannot exceed the observed 
value of 5T/T ~ 5 x 10" 6 03]. 

A detailed analysis of the CMB bound yields [H EH 

Q g h 2 10 o < 7 x lCT 11 (^) (H Q < f < 30 x H ) (3.42) 
where Q g = fP(f)/p , p = 3cm 2 pl Hl / '(8irh) and H = h 100 x 100A;m/(s x 
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Mpc) with h 100 = 0.7. The CMB bound for the spectrum ()3.34j) evaluated 
at 2-kHq Hertz reads 2 

/ \ 41+4 / \ 2 / jj \ 3 / rr \ 41-1 

l>(27r)- 4 '/ 2 - 4 '(2/ 2 -3/+l) 2 (^) ( ^r-r) — 



(3.43) 



and consequently 



f(l, H 1} r) = -107.69+/ ( 28.10 + 2 log 10 ( J J -(2Z-2) log 10 ( 1 + 



y y v« ^ ■ j 

(3.44) 

+ (-4/ + 2) log 10 I + 2 log 10 (2/ 2 - 3/ + 1) < 0. 

We next consider different values for I and r. 
(i) When I = 1.1, the relation (|3~4ljl reads 

/(l.l, F l5 r) = -76.71 + 2.20 log 10 (J^j - 0.2 log 10 (l + ^p#i"r) < 0, 

(3.45) 

see Fig. 13.31 From it we observe that: 

1. For r < r^ 1 " 1 = 1.22 x 10~ 35 s, the condition ((S3HD is satisfied if #i < 
8.13 x 10 34 s _1 and conflicts with ()3.38jl . which in the most favorable 
case is Hi = 8.13 x 10 34 s _1 for r = rjr 1 - 1 . For I = 1.1, there is 
no compatibility with the observed CMB anisotropy when r < rjr 1 ' 1 . 
Thus, this range of r is ruled out. 

2. For r > rjr 1 ' 1 , one obtains -f/i < H 1 c =1a (t) from the condition (|3.45j) . 
H 1 c =1a (t) is always larger than r -1 in the range considered, e.g. H 1 c =1a (t 
10~ 30 s) = 2.40 x 10 35 s _1 . Taking into account (|3~38|) one obtains that 
the condition (J3.45j) is satisfied for r _1 < Hi < _£^ =l1 (t). 



(ii) When I = 2 (the extreme case in which the expansion is entirely 
dominated by the MBHs) we have 

f{2,H h r) = -51.89+4 log 10 (J^J -21og 10 (l + < 0, (3.46) 

see figure l3~B1 Inspection of ()3.46|) and Fig. 13.31 reveals that: 

2 It is necessary to multiply H3.34|l for h/c 3 in order to obtain the right units. 
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a b c a b c 




log .= CH i' E) i°g„,(H,s} 



Figure 3.3: The left panel depicts /(l.l, r)t>s. log 10 ifi for a) r = 10 40 s, 
b) r = r^ 1 - 1 = 1.23 x 10~ 35 s, and c) r — 10~ 30 s. The right panel depicts 
f(2, H x , t) for a) r = l(T 20 s, b) r = r^f 2 = 6.75 x l(T 14 s, and c) r = l(T 10 s. 
Conditions (|3.44j) and ()3.38|) are satisfied for certain ranges of r and H\ in 
each case. 

1. For r < t 1 = 2 = 6.75 x 10" 14 s, one obtains Hx < 10 13 s 1 which is totally 
incompatible with condition (|3.38)L Hx > 1.48 x 10 13 s _1 for r = r l c =2 
in the most favorable case. Thus, the region r < t 1 c =2 is ruled out as 
predicts an excess of anisotropy in the CMB. 

2. For r > r^ 2 , one obtains Hx < H 1 c =2 (t) from the condition (J3.46|) . 
H 1 c =2 (t) is always larger than 1.48 x lO 13 ^" 1 for r in the range consid- 
ered, e.g. H 1 = 2 (t = l(T 10 s) = 1.35 x lO 16 ^ 1 . Conditions (EPSl) and 
()3.46|) are both satisfied in this range for r _1 < H\ < H 1 c =2 (t). 



We may conclude by saying that the condition ()3.44|) leads to different 
allowed ranges for Hx and r for each I considered, although their interpreta- 
tion is rather similar. For r < t 1 c the condition of minimum duration of the 
"MBHs+rad" era ()3.38|) and the upper bound given by the CMB anisotropy 
are incompatible. However, for r > t 1 c these two conditions are compatible 
for r- 1 <Hx< H l c (r). 
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3.3 Conclusions 

We have calculated the power spectrum of GWs in a universe that begins 
with an inflationary phase, followed by a phase dominated by a mixture of 
MBHs and radiation, then a radiation dominated phase (after the MBHs 
evaporated), and finally a dust dominated phase. The spectrum depends 
just on three free parameters, namely H\ the Hubble factor at the tran- 
sition inflation- "MBHs+rad" era, r, the cosmological time span of the 
"MBHs+rad" era, and the power /, being a(r/) oc rf the scale factor of the 
"MBHs+rad" era with 1 < I < 2. 

The upper bound on the spectrum of GWs obtained from the CMB 
anisotropy places severe constraints on Hi and r. For each value of I consid- 
ered, there is a minimum value of r, t 1 c , compatible with the CMB anisotropy. 
There is a range of r, r > r l c , for which r _1 < Hi < H 1 c (t) satisfies the CMB 
upper bound. 

The four-stage scenario predicts a much lower power spectrum of GWs 
than the conventional three-stage scenario. If LISA fails to detect the GWs 
spectrum at the level predicted by the three-stage scenario, a possible ex- 
planation might be that the "MBHs+rad" era took place. Likewise, once 
the spectrum is successfully measured we will be able to learn from it the 
proportion of MBHs and radiation in the mixture phase as well as the time 
span of this era. 



Chapter 4 



Gravitational waves and 
present accelerated expansion 

In this chapter we show how the power spectrum as well as the dimensionless 
density parameter of the GWs created from the initial vacuum state may help 
present (and future) observers ascertain whether the expansion phase they 
are living in is accelerated or not and if accelerated, which law follows the 
scale factor [45] . The latter would facilitate enormously to discriminate the 
nature of dark energy between a large variety of proposed models (cosmo- 
logical constant, quintessence fields, interacting quintessence, tachyon fields, 
Chaplygin gas, etc) jj^]- To this end we calculate the power spectrum and 
energy density of the GWs when the transitions to the dark energy era and 
second dust era are considered. Obviously, the latter power spectrum lies at 
the future and depending on the model under consideration it may take very 
long for the Universe to enter the second dust era. 

4.1 Accelerated expansion and decaying dark 
energy 

Nowadays the observational data regarding the apparent luminosity of su- 
pernovae type la, together with the discovery of CMB angular temperature 
fluctuations on degree scales and measurements of the power spectrum of 
galaxy clustering, strongly suggests that the Universe is nearly flat and that 
its expansion is accelerating at present J3HJ HI] • In actual fact the debate 
now focuses on when the acceleration did really commence, if it is just a 
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transient phenomenon or it is to last forever, and above all the nature of the 
dark energy. 

In Einstein gravity the accelerated expansion is commonly associated to a 
sufficiently high negative pressure which might be provided by a cosmological 
constant A (vacuum energy), whose equation of state is p\ = — p\, with 
Pa — A/(8?r) = constant. The observational data seems to indicate that 
the cosmological constant contributes about the 70% of the energy of the 
Universe, meanwhile the remaining 30% comes from non-relativistic matter 
(i.e., dust). Thus, the question arises: "Why are the vacuum and matter 
energy densities of precisely the same order today?" , which is known in the 
literature as the coincidence problem |IH] . 

A possible answer to this question considers that the acceleration of the 
Universe is associated to a sort of dynamical energy, the so-called dark en- 
ergy, that violates the strong energy condition and clusters only at the largest 
accessible scales In such the present state of the Universe would 

be dominated by dark energy and since it redshifts more slowly with expan- 
sion than dust, the contribution of the latter is bound to become negligible 
at late times. 

In an attempt to evade the particle horizon problem posed by an everlast- 
ing accelerated expansion to string/M type theories [49], some models pro- 
pose dark energy potentials such that the current accelerated phase would 
be just transitory and sooner or later the expansion would revert to the 
Einstein-De Sitter law a(t) oc t 2 ^ 3 , thereby slowing down (second dust era) 
|50j . The possibility that dark energy could be unstable is in fact suggested 
by the remarkable qualitative analogy between the presence of dark energy 
today and the properties of a different type of dark energy - the inflaton field 
- postulated in the inflationary scenario of the early Universe -see e.g., |5T] . 
This analogy has two main points. On one hand, it makes natural that a 
form of matter with negative pressure could have dominated the Universe in 
a distant past, since a similar form of matter dominates the Universe today. 
On the other hand, as the dark energy in the early Universe was unstable 
and decayed aeons ago, one might be tempted to ask whether the nature of 
dark energy observed today would be any different. 

For our purposes, we shall assume a simplified model of decaying dark 
energy in which the usual dust dominated stage is followed by an accelerated 
era where the adiabatic index 7 of the fluid that dominates that era is a 
constant that lies in the range [0,2/3). Posteriorly, the dark energy decays 
in a time span much lower than the duration of the accelerated era and 
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the Universe resumes the decelerated expansion dominated by the cold dark 
matter. Models with dark energy whose evolution mimics that of dust can be 
also taken in consideration in our description as they are formally equivalent 
to the model of above. 



4.2 Bogoliubov coefficients in the decaying 
dark energy scenario 

In this section we evaluate the coefficients of Bogoliubov in the simplified 
model previously suggested, a spatially flat FRW scenario initially De Sitter, 
then dominated by radiation, followed by a dust dominated era, an accel- 
erated expansion era dominated by dark energy, and finally a second dust 
era. 

The scale factor in terms of the conformal time reads 

i- (— oo < 77 < rji < 0), De Sitter era 



a(rj) 



H 1V 

77 — 7 (r] — 2t]i) (r/i < r/ < 172), radiation era 



Ilk? {V+ m-Z lY (m<r]<m), first dust era 

* (??3+??2-4?7i) 2 -' 
AH^ri[(r)2-2r)x) 

1 7T-777 



1 



(5) ^t^ffe-^i) fa)' i.V3<V<Vd, dark energy era 

?f ((I4H4) 2 [r] — jj- ) (774 < r/), second dust era 



(4.1) 

where I < —1, r\i — 77 + | [(—2// + 1) r/ 3 + r/ 2 — 4^], the subindexes 1, 2, 3, 4 
correspond to sudden transitions from inflation to radiation era, from radi- 
ation to first dust era, from first dust era to dark energy era and from the 
latter to the second dust era, respectively; Hi is the Hubble factor at the 
instant r\ = r/i. The present time rj lies in the range [773,774], it is to say in 
the dark energy dominated era. 

GWs which are inside the Hubble radius have a wave number k lower than 
a(r])H(i]). Figure l4~T1 sketches the evolution of a(i])H(r]) in this scenario. 
During the inflationary and dark energy eras a(i])H(r]) increases with 77, and 
decreases in the other eras. As a consequence, 04^4 results higher than 03-^3. 
Choosing /, 773 and 774 in such a way that 

\ -i/' / \ 1/2 

O4 \ / 02 



- >1, (4.2) 

0-3/ V°3/ 
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we have that a^H^ is also higher than a 2 H 2 . We assume that a^H^ is lower 
than a\H\ throughout. 

The scale factor (j4.1j) formally coincides with these of the standard three- 
stage model given by (|3.1|) until the instant r] 3 . Thus, the Bogoliubov coeffi- 
cients «i and Pi are equal to those of equation (|3.5|) in the range k < 2-aa\Hx 
and they are 1 and for k > a\H\, respectively. Analogously, a 2 and [3 2 
coincides with those of ()3.8|) for k < 2na 2 H 2 and are 1 and for k > a 2 H 2 . 

In the dark energy era the solution to Lifshitz's equation (j2.10J) reads 

IH = (V^me^k^xy'H^ixt), (4.3) 
where x\ = krji. The Bogoliubov transformation 

Ml) = okM + Psrtiv) (4-4) 

relates the dust and dark energy modes. 

Using the well-known relation for Hankel functions jHHl 

H\ 2 \ /2 {x) ~ e " lX \/^^T(-2)^(l + «), (4-5) 
valid when i< 1 and I < — 1, it follows that 

as = (6(-l)- / /^)^- 2 + ^(^ 1 ), 
/3a = -(6(-l)-¥^)^- 2 + ^- X ), 

for k < 2ira(rj 3 )H 3 and 013 = 1, ^3 = for k > 2Tra(r] 3 )H 3 , where H 3 is the 
Hubble function evaluated at r/ 3 and rj m3 = r/ 3 + r\ 2 — 4r/i. The transition 
between the first dust era to the dark energy era has a time span of the 
order H 3 X . This is much shorter than the period of the waves we are consid- 
ering and therefore it may be assumed instantaneous when calculating the 
coefficients. As a consequence, the time span from this transition till today, 
r = to — t 3 , must be larger than H 3 X , otherwise the transition first dust 
era-dark energy era would be too close to the present time for our formalism 
to apply. This condition places a lower bound on the value of the redshift — . 
When I = -1 the condition is ^ > 2.72, ^ > 2.25 when I = -2, ^ > 2.15 
when I = —3, and — > 2 when / — > —00. This bound is compatible with 
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Figure 4.1: Evolution of a{rj)H(rj) in a universe with scale factor given by 
Eq. (j4.1|) under the assumption a±H± > a 2 H 2 . The quantity r]f is defined 
as the instant in the dark energy era in which 0(77/0) #(77/0) = 02-^2, Vfi as 
the instant of the second dust era such that a(?7/i)if (77/1) = 02-^2, and 77/2 
as the instant of the second dust era in which 0(77/0)^(77/0) = CI3H3. 
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the accepted values for ^ £ [1.5, 11] (see e.g., |S21). Henceforward we will 
consider ^ larger than 2.72 and no larger than 11. 

The solution to Lifshitz's equation in the second dust era (i.e., the one 
following the dark energy era) is 

fi SD = {^/2)e^k~ 1 / 2 xf D Hfl{x SD ), (4.6) 

where x SD = k (rj - 7/4 + ^ J . 

The Bogoliubov coefficients relating the modes of the dark energy era 
with the modes of the second dust era are 

M = {v) + Piv*sD (V) • (4-7) 

Likewise, the continuity of fi at 774 implies 

«i = (-l) l l2 l t^x\f[l 2 + il(l-2)x H + 0(xl)], 

A = (-l)^2'^4- 2 [/ 2 + z/(/ + 2)x i4 + 0(4)], 

for k < 2na i H A) and a\ — 1, /3| = for A; > 2na i H i) where i?4 is the Hubble 
function evaluated at the transition time 774. 

4.3 Power spectrums 

In this section we evaluate the total coefficients and calculate the current 
power spectrum as well as the power spectrum in the second dust era. We 
find that the current power spectrum coincides with those of the three-stage 
model of section 13.11 but it evolves differently. We also find two possible 
shapes for the power spectrum in the second dust era, whether condition 
flOJ) is fulfilled or not. 

4.3.1 Current power spectrum 

To evaluate the present power spectrum one must bear in mind that a^H^ > 
(I3H3 which implies the wave lengths of the perturbations created at the 
transition dust era-dark energy era are larger than the present Hubble radius 
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H 1 [SHj- One must also consider the possibility that aoH > o, 2 H 2 , it is to 
say 



For I 



-1 and assuming — 



a 3 



(4- 



10 4 [SB] this condition implies ^ > 21.5 



^ > 100 when I = -2, and a > 251.2 when Z = -3. The values for 2a 

as ' a3 a3 

considered by us are larger than 2.72 and no lower than 11, consequently we 



can safely assume a H < a 2 H 2 . 

For k > 2ira\Hi, we find that a>Tr2 ■ 
k > 2na 2 H 2 , the coefficients are ax r 2 
k < 2Tra(r] 2 )H 2 we obtain j271l2^j 



1, Ptt2 = 0; in the range 2-na\H\ > 
- a 1 and Ptt2 = Pi, and finally for 



P 



'Tr2 



8k 3 r]fr] R2 ' 



(4.9) 



Thus, the number of GWs at the present time rj Q created from the initial 
vacuum state is (N^) = \PTr2f ~ t ^~ 6 ( ? 7o) f° r ^(^0) < 27r(a 2 /ao)i?2, a;_4 ( ? 7o) 
for 2Tr(ai/ao)Hi > u(r]o) > 2ir(a 2 /ao)H 2 , and zero for u(r] ) > 2^{a\/ a^Hi, 
where we have used the present value of the frequency, ^(770) = k/a . 

In summary, the current power spectrum of GWs in this scenario is 



P(uj) 







(a? (770) > 2Ti(a 1 /a )Hi) 



1 

4tt 2 



(2ir(a 2 /a )H 2 < uj(r) ) < 2ir(ai/ao)H 1 ) 



16vr 2 



-3 



(2irH < uj(r}o) < 2n(a 2 /a )H 2 ) . 

(4.10) 

While this power spectrum is not at variance with the power spectrum of 
the conventional three-stage scenario evaluated at section I3~T1 it evolves dif- 
ferently. The power spectrum in the dark energy scenario at T] — 773 formally 
coincides with Eq. (|4.10|) but with a 3 substituted by ao throughout, and 
from then up to now waves with 2irazH^ < k < 27rao-£/o cease to contribute 
to the spectrum as soon as their wave length exceeds the Hubble radius. By 
contrast, in the three-stage scenario waves are continuously being added to 
the spectrum. As we shall see in section 14 A\ this implies that the evolution 
of the energy density of the gravitational waves in the three-stage scenario 
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differs from the scenario in which the Universe expansion is dominated by 
dark energy. 



4.3.2 Power spectrum in the second dust era 

Here we evaluate the power spectrum at some future time rj larger than 
77/2 for which the waves created at the transition dust era-dark energy era 
(V — 773) are considered in the spectrum by the first time (see Figure l4~TJ) . 
Let «Tr4 and (3t T 4 be the Bogoliubov coefficients relating the modes of the 
inflationary era to the modes of the second dust era. Because of condition 
(j4.2j) we must consider two possibilities with two different power spectrums. 

(i) If condition ()4.2j) is not fulfilled, then aiH x > a 2 H 2 > a^H^ > a 3 H 3 
and the power spectrum can be obtained from the following total coefficients. 
In the range k > 2tt a\H\ the total coefficients are «Tr4 = 1 and Ptta = 0. For 
2naiH 1 > k > 2na 2 H 2 , the coefficients are «T r4 = a% and /3t>4 = Pi, where 
ai and (3\ are defined in Eq. ([3.5)1 . For 2-aa 2 H 2 > k > 2-aatJi^ oiTri = &Tri 
and /?Tr4 = Ptt2 where /3t T 2 is defined in Eq. (|4.9|) and the dominant term 
of aTr2 coincides with Ptt2- 

When 2'Ka^H^ > k > 2na 3 H 3 , except for a 3 and /3 3 , all the coefficients 
obtained in the previous section must be considered when evaluating «y r4 
and Ptt±, therefore 



a Tr4: — 0(\oiT r 2 + P\ Ptt2 



PttA ~ Pi a Tr2 + °4 Ptt2 



32 1 j (-/)! 



A;'- 5 (4.11) 



rfim2ml +2 



32 1 ' (-1)1 J vfvmvr^ 2 



k l - 5 (AA2) 



Finally, for 2na 3 H 3 > k > 27rH(r]), we get 



P\ 



Tri 



16 



-l l+2 2 l 



-2l)\ 



21 



,.21-5 



vfVR2vJ 3 Vu +2 



(4.13) 
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( (w{rj) > 2ir(a 1 /a(ri))H 1 ), 

4 



P(u) 



(^y) Hfu- 1 (2n(a 2 /a( V ))H 2 < u(r,) < 27r(a 1 /a(r 7 )) J ff 1 ) , 



16-7T 2 V a2 



{2^(0,4/ a{rj))H 4 < w(t/) < 2Tr(a 2 /a(r]))H 2 ) 



< it 1 ) 



2/-10 



16-4/ 



/-4 



1-4+4/1 



-8+21-4/1 



Hl°~V- 7 



(2w(a 3 /a(r)))H 3 < u){rj) < 27r(a 4 /a(r?))iJ 4 ) 



27^4Z 2 4«-6 / (~2Z)! \ 4 / ai 
7? V ( - ! / \ a ° 



a (v) \ 

a± J 



4Z — 10 



16-8/ 



2Z-4 



2Z-4+4/Z 



/ fly 



41-8-4/1 



10-4/, ,4/-7 



(2irH( V ) < 0,(77) < 2tt {a 3 /a{rj))H 3 ) . 

(4.14) 

(ii) If condition ()4.2)1 is fulfilled, then aiiJi > a 4 if 4 > a 2 H 2 > a 3 H 3 . As in 
the previous case, in the range k > 2-na\H\ the total coefficients are «Tr4 = 1 
and /?ty4 = 0. For 2-Ka\H\ > k > 27ra 4 if 4 , the coefficients are ar r 4 = «i and 
Ptt4 — Pi- In the range 27ra 4 if 4 > k > 2na 2 H 2 , we obtain 



<4>4 = a[a, + ft % ~ ^ (-1)-' M \ 2 / )! - 1 m k 1 ^ 



Ptt4 = P 4 ai + a 4 Pi / 2 2 . 

o I 'J- %?7j 4 



Again, for 2ira 2 H 2 > k > 2ira 3 H 3 , the total coefficients are given by Eq. 
(j4~TT]) while for 2na 3 H 3 > k > 2nH(r]) they obey Eq. (jO^jl . The power 
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spectrum in this case is 

f (u(ri) > 27r(a 1 /afa))ffi) , 
l 



P{w) ~ < 



^(^y) fffr" 1 (27r(a4/o(f7))ff4 < w(t/) < 27r(ai/o(f7))ffi) 



Itt 



9 /22+2022-6 / (-22)! \ / oi 



a 4 I 



22-6 



-12 



2-1 



1-3+2/1 



21-4-4/1 



6-21, ,21-3 



(2TT(a 2 /a(r]))H 2 < u(rj) < 2ir(ci4/a(r]))H i 



9 /22r>22-8 ( (— gOj 

^ 2 l-R)T 



a.4 / 



27^4« 2 4«-6 ( (-2Q! 



22-10 



16-42 



ai 
ao 

10-22, ,22-7 



2-4 



2-4+4/2 



22-8-4/2 



(2tt (a 3 /a(rj))H 3 < uj{r,) < 2Tr(a 2 /a(ri))H 2 ) , 



FIT 

42-iO 



16-82 



22-4 



22-4+4/2 



42-8-4/2 



#10-42^42-7 



(27tH( V ) < u(rj) < 2ir{a 3 /a{ri))H 3 ) . 

(4.15) 

The power spectrum governed by Eq. (|4.15jl is plotted in Fig. 14 .21 for different 
choices of the free parameters 1, —, — and as well as the power spectrum 
assuming the three-stage model, i.e., non-accelerated phase and no second 
dust era. 

The shape of the power spectrum given by (|4.15jl in the range 2tcH(7]) < 
oo (rj) < 27r(as/a(r]))H 3 is the same in both cases as in this range all the 
coefficients are present in the evaluation of the total coefficients. It is inter- 
esting to see how markedly this spectrum differs from the one arising in the 
three-stage model (dot-dashed line) at low frequencies. 



Topological defects 

Up to now we assumed that at the end of the dark energy era the Universe 
will evolve as if it became dominated by dust once again. Nevertheless if 
the expansion achieved in the accelerated phase were large enough, either 
cosmic strings, or domain walls, or a cosmological constant will take over 
instead. We will not consider, however, cosmic strings (whose equation of 
state is p% = —\pk) for, as pointed out by Maia [21], it seems problematic 
to define an adiabatic vacuum in an era dominated by these topological 
defects since the creation-annihilation operators, A^, Atx, fail to satisfy 
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Figure 4.2: Power spectrum given by Eq. ()4.15|) for different values of I. The 
left panel depicts the power spectrum at the instant rj for which = 10 22 
(when fig (77) is near unity), with ^ — 11, and ^ = 10 6 (to fulfill condition 
(|4.2j) ). The right panel shows the power spectrum at the instant rj for which 
^ = 10 16 (again when tt g (r]) is near unity), with 2a = n ; an d M = 10 6 . 
For the sake of comparison the power spectrum of the three-stage scenario 
(De Sitter inflation, radiation, dust with / = 2) is also shown in both panels 
(dot-dashed line). Notice the difference in slopes at the dust era. 
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the commutation relations (i.e., condition (J2.29j) ) in the range of frequencies 
where one should expect GWs amplification. In short, our approach, as it is, 
does not apply to this case. 

As for domain walls (topological stable defects of second order with equa- 
tion of state pd w = —\pdw and energy density that varies as a~ 1 (^) -see e.g., 
ED OH] ) j once the dark energy evolved as pressureless matter at r\ = 774 the 
scale factor may be approximated by 



a(rj > 774) = 4a 4 (a 4 H 4 ) yV + Vi + ^jfj > ( 4 - 16 ) 

so long as a -1 ^> a~ 3 . That is to say, for 77 > 7/4 the expansion of the Universe 
is again accelerated whereby a(rj)H{rj) resumes growing. The GWs will be 
leaving the Hubble radius as soon as H~ l becomes smaller than their wave 
length, and eventually none of them will contribute to the spectrum. 

Finally, we consider the existence of a positive cosmological constant A. 
(Recall that p\ = A/ (8tt) and p\ = —Pa). Once the dark energy dynam- 
ically mimicked dust the Universe will become dominated by a very tiny 
cosmological constant. The corresponding scale factor is 

/ 1 \ - 1 

r-l 



a(rj > t/ 4 ) = # 4 ( -77 + V* + ^jfj ■ (4-17) 

Once again, the expansion is accelerated and GWs will leave the Hubble 
radius and, in the long run, none of them will contribute to the spectrum. 



4.4 Energy density of the gravitational waves 



Now we are in position to evaluate the energy density of the GWs in terms 
of the conformal time by integrating the power spectrum P(u) obtained in 
the previous section. As we shall see, the evolution of the energy density 
strongly depends on the free parameters of the model. 

Its current value, evaluated from Eq. (j4.10j) . can be approximated by j2Z] 

P 9 (Vo) - ^ (^X Ht{2*H Q y\ (4.18) 

d,2-K l c 6 \a 2 J \a J 

where, in this section, we return to conventional units. 
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To study the evolution of p g {rj) from this point onward we first consider 
that a 4 if 4 > a 2 if 2 - In this case, p g {jf) evolves as (J4.18|) with H{rf) and a{rj) 
substituted by H and ao, respectively, till some instant 77/0 in the range 
Vo < Vfo < When 77 ^ 77/0 the GWs with u>(r]fo) < 2n (02/0(77/0)) H 2 
must no longer be considered in evaluating p g as their wave length exceed 
the Hubble radius. Consequently 

h ( ai \ 4 T 4 / diHx 



P 9 (Vfo < V < m) - ^ Hfln ) . (4.19) 

For 77 = 774, the gravitational waves created at the transition dark energy 
era-second dust era begin contributing to p g thereby, 



a o / \ a o J \ a o ) V a 4 / — 2/ + 2 

where p g {j]^) corresponds to Eq. (j4.19|) evaluated at 77 = 774. For 77 = 77/1 > 
774 where 77 is defined by the condition a(r)fi)H(r]fi) = a^H^ (see Figure 
I4.ip . the gravitational waves which left the Hubble radius at 7?/ reenter it, 
therefore 

< ^^) = ^)(^) 4 + ^ a - 8 (^) 2 © 16 " 1 '(^) 

>V- 4 /a3Y- 1+4/ ' /M 2 '- 8 -" /«Wf^ 1, i[1M (2>rff(-,))« 



oo/ V a o/ \aoJ V a 4 / 1 -2/ + 6 

where p g {r]fi) corresponds to Eq. (|4.20|) evaluated at 77 = 77/1. For 77 = 77/2, 
with 77/2 defined by a (77/2) -£^(77/2) = 03-^3 (see Figure l4~T|) . the gravitational 
waves created at the transition first dust era-dark energy era have wave 
lengths shorter than the Hubble radius for first time, and from that point on 
the density of gravitational waves can be approximated by 

«/2 V , 27» l4Jo4J _ 6 /(-20!\V«i xl,; "' S/ 



P 9 (V > m )*p, M ^) + ^«-^J (4.22) 



a o J \ a o J \ a o J V a 4 / — 4/ + 6 



where p g (/7/ 2 ) corresponds to Eq. (|4.21j l evaluated at 77 = 77/2. 
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In the case that CI4H4 > (I2H2, Eq. 1)4.18)1 dictates the evolution of p g 
between 770 till 774. Then, from r/ 4 till 77/2, Pg obeys Eq. ()4.21j) (note that 77/1 
cannot be defined in this case) and from 77/2 onwards p g obeys Eq. ()4.22|) . 

A natural restriction on p g (f]) is that it must be lower than the total 
energy density of the flat FRW universe 

P(V) = j^^U (4.23) 



where rapi — ^hc/G stands for the Planck mass. 

It seems reasonable to consider Hi ~ 10 38 s _1 as it corresponds to the 
grand unification energy scale of inflationary models 0H1 E3] . The redshift 
^, relating the present value of the scale factor with the scale factor at the 
transition radiation era-first dust era, may be taken as 10 4 |38j . and the 
current value of the Hubble factor H is estimated to be 2.24 x lCT 18 ,?" 1 
[HI EE]- It then follows 



^ = 1.50 xlW^ , H 2 = 2.24 x 10- 12 ( * ) 

a \a 3 J \a 3/ 



2 1 1 



1+4 / \ -1 



1 



H, II [-) Y\ //, //o(-) 

a 3 J \a J 

H„ = 2.24 x loW (^V' + V, 



, 3/2 

04 



! /2 

where we have used the relation 

"^ri'KW)) H " (424) 

valid in the second dust era (see Eq. (|4.1)l ) with afiHfi = (I2H2 and af 2 Hf 2 = 
CI3H3 to evaluate Hfi (only defined if a^H^ > 02-^2) and if/2, respectively. In 
our model, there are only three parameters, namely I < —1, 2.72 < — < 11, 
and ^. 

ao 

We are now in position to evaluate the evolution of the dimensionless 
density parameter £l g (rj) = p g (r]) / p(r)) . Its current value is [2JJ 

1 Gh f ciq \ 2 f (ii x 8 
48tt 3 c 5 \a 2 J Vao 



^ - 7^-^ ( - ) ( - ) H?H^\ (4.25) 
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which in our description boils down to 



n g0 ~ 2.00 x 10" 13 I — 



a 3 



(4.26) 



flgo is much lower than unity for any choice of ^ and I in the above intervals 



L g0 

At later times Q g evolves as 



ttgiv > Vo) = ^ 



90 



0(77) 
a 



(4.27) 



Hn. It is obvious that 



where we have used H(r]4 > rj > rjo) = (j^- 
Q g {r]) is a decreasing function of 77. For rj > r] Q we shall distinguish the two 
cases mentioned in the previous section. 

When condition (|4.2j) is fulfilled, the evolution of fl g is given by Eq. ()4.27|) 
till r] = rjfQ. Then, p g changes in shape from 77 — rjfo till 77 = 774, as we have 
seen. Consequently 



fi 9 (?74 > V > Vfo) 



2 Gh ( ao 



3ir c 5 \a 3 



a 2 / 0(77) 



(4.28) 



x In 



03 ^ 

ao. 



1 _ j_ 

4 21 



a 



is a decreasing function of 77. Finally, from 774 on, Q g evolves as 



04 
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mv)) 



21-4 



Vg(Vf2 > 77 > 77/1) ~ n g (77/1) 
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ao/ 
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As follows from ()4.24j) and (J4.29j) - (j4.31j) . in each case, the first term redshifts 
with expansion while the second term of Q g grows with rj (recall that I < — 1). 
Conditions Q g (r]fi > i] > 7/4) <C 1 and fl g (rjf2 > r\ > i]fi) <C 1 become just 
fl g (rj = rjfi) ^ 1 and £l g (rj = 77/2) <C 1, which are true in either case whatever 
the choice of parameters. Finally, from Eq. (J4.31|) we may conclude that at 
some future time larger than 7iy 2 the condition fi 9 < 1 will no longer be 
fulfilled and the linear approximation in which our approach rests will cease 
to apply. 

In the opposite case, when condition ()4.2|) is not fulfilled, fl g also grows 
following Eq. ()4.27j) till 77 = 7/4, in the interval 7/4 < rj < 77/2 it grows 
according to Eq. ()4.30j) . and according to Eq. (j4.31|) from 7// 2 onwards. Our 
conclusions of the first case regarding the evolution of Q g during these time 
intervals still hold true. 

The behavior of the density parameter Q g differs from one scenario to 
the other. In the three-stage model (inflation, radiation, dust) Q g remains 
constant during the dust era [2ZJ. I 11 a four-stage model, with a dark en- 
ergy era right after the conventional dust era, fl g sharply decreases in a I 
dependent way during the fourth (accelerated) era since long wave lengths 
are continuously exiting the Hubble sphere jHSj, see Fig. 14.31 

By contrast, if a second dust era followed the accelerated (dark energy) 
era, fl g would grow in this second dust era because long wave lengths would 
continuously be entering that sphere, see Fig. 14.41 This immediately suggests 
a criterion to be used by future observers to ascertain whether the era they are 
living in is still our accelerated, dark energy-dominated, era or a subsequent 
non-accelerated era. By measuring Q g at conveniently spaced instants they 
shall be able to tell. Further, if that era were the accelerated one and the Q g 
measurements were accurate enough they will be able to find out the value 
of the parameter I occurring in the expansion law given by Eq. (j4.1|) . The 
lower I, the higher the slope of Q g in the second dust era. 

One may argue, however, that these observers may know more easily 
from supernova data. Nonetheless, if this epoch lies in the faraway future it 
may well happen that by then the ability of galaxies to generate stars (and 
hence enough supernovae) has seriously gone down and as a result this prime 
method might be unavailable or severely impaired. At any rate, even if there 
were plenty of supernova, the simple gravitational wave method just outlined 
could still play a complementary role. 
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Figure 4.3: Left panel: Evolution of the density parameter Q g with the scale 
factor from the beginning of the radiation era till the present time. The 
solid line shows the density parameter predicted by the three-stage model 
(De Sitter inflation, radiation, dust). The dashed and dotted lines depict 
the density parameter predicted by a four-stage model having a dark energy 
era right after the dust era with / = —3, —1, respectively, and a /a 3 = 11. 
Note that the four-stage model predicts a lower present value for f2 9 than the 
three-stage one as fl g decreases, in a I dependent fashion, during the dark 
energy era. The right panel is a blow up of the region in which the four-stage 
models (/ = —3, —1), notably differ from the three-stage one. 
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Figure 4.4: Evolution of the density parameter Q g given by Eqs. (|4.2fij) - (|4.31jl 
with the scale factor in a five-stage scenario (De Sitter inflation-radiation- 
dust-dark energy-second dust era) from the beginning of the dark energy 
era. The second dust era starts at the time in which the graphs attain their 
minimum. The solid, dotted and dashed lines correspond to I = —1, I = —1.5 
and / = —2, respectively. The density parameter decreases in the dark energy 
era and increases in the second dust era till its value is comparable to unity. 
From that point on our approach (being linear) ceases to apply. We have 
chosen a dark energy era which begins when a 3 /a = 1/11 and ends when 
a 4 /a = 10 6 . 
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4.5 Conclusions 

In this chapter we studied the power spectrum and the energy density evolu- 
tion of the relic gravitational waves generated at the big bang by considering 
the transitions between successive stages of the Universe expansion. In par- 
ticular, we considered the effect of the present phase of accelerated expansion 
as well as a hypothetical second dust phase that may come right after the 
present one. As it turns out, the power spectrum at the current accelerated 
era, Eq. (j4.1Uj) . formally coincides with the power spectrum of the conven- 
tional three-stage scenario. As a consequence, measurements of P{ui) will 
not directly tell us if the Universe expansion is currently accelerated (as we 
know from the high redshift supernove data) or non-accelerated. However, 
the density parameter of the gravitational waves evolves differently during 
these two phases: it stays constant in the decelerated one and goes down in 
the accelerated era in a I dependent manner. Therefore, the present value of 
Q g may not only confirm us the current acceleration but also may help de- 
termine the value the parameter I -see Fig. 14.31 - and hence give us invaluable 
information about the nature of dark energy. 

In the faraway future measurements of P{uj), if sufficiently accurate, will 
be able to tell if the Universe is still under accelerated expansion (driven 
by dark energy) or has entered a hypothetical decelerated phase (second 
dust era) suggested by different authors jSU]. This may also be ascertained 
by measuring the density parameter of the gravitational waves at different 
instants to see whether it decreases or increases. 
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Chapter 5 



Gravitational waves entropy 
and the generalized second law 

In this chapter we study the evolution of the entropy associated to the GWs 
as well as the generalized second law (GSL) of gravitational thermodynamics 
in the present era of accelerated expansion of the Universe. In spite of the 
fact that the entropy of matter and relic gravitational waves inside the event 
horizon diminish, the mentioned law is fulfilled provided the expression for 
the entropy density of the gravitational waves satisfies a certain condition 
|57j. Section I5~D gives the power spectrum of the GWs at the beginning 
of the present era of accelerated expansion. In section 15.21 the entropy of 
the GWs is defined and its evolution during that era is found. Finally, in 
section 15. HI an upper bound on the entropy density of the GWs is obtained 
by straightforward application of the GSL. 

5.1 Gravitational waves in the dark energy 
era 

As explained in the previous chapter, the current era of cosmic acceleration 
is believed to be dominated by the dark energy. For our purposes in this 
chapter, we shall assume the Universe is currently dominated by a form of 
everlasting dark energy of constant 7 lying in the range (0, 2/3) (i.e., phantom 
energy and cosmological constant dominated universes are excluded). The 
dependence of the scale factor on the conformal time is formally equal to the 
scale factor in (|4.1jl . without the second dust era, i.e., 
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(iji < v < m), 

(r] 2 < r] < rj 3 ), 
(V3 < V), 



De Sitter era 
radiation era 

dust era 
dark energy era 
(5.1) 



where I < — 1, 



, Vi — V + | V3 + V2 — 4r]i] . Again, the subindexes 1, 2, 3 

correspond to sudden transitions from De Sitter era to radiation era, from 
radiation to dust era, and from dust era to dark energy era, respectively, and 
Hi is the Hubble factor at the instant rj = rji. The Hubble factor during the 
current dark energy era obeys 



H(rj) 



«3 

a{rj) 



(5.2) 



The evolution of the quantity aH in terms of the conformal time rj is sketched 
in Fig. 15.11 

The modes solution to the gravitational wave equation during the De 
Sitter era are related with those of the final dark energy era by a Bogoliubov 
transformation with coefficients ax r 2 and /3ty2, which are formally equal to 
those found in subsection 14.3.11 At the beginning of the dark energy era, 
rj = 773, the power spectrum was 



' (w(7fe) > 27r(ai/a 3 )i?i) , 

P(u) ~ \ h {f)' H t^ 1 (2n(a 2 /a 3 )H 2 < uj(rj 3 ) < 2n(a 1 /a 3 )H 1 ) , 

16^ ( j) 2 (S) 8 H ^ < u fa) < Ma2/a 3 )H 2 ) . 

(5.3) 

During the radiation and dust eras a(rj)H(rj) decreases with rj and in- 
creases during the De Sitter and dark energy eras. Consequently, as we 
explained before, GWs are continuously leaving the Hubble radius during 
the accelerated dark energy era ES]- At some instant 77/0, defined by 
a(r)fo)H(r)fo) = a 2 H 2 , the third term in (j5.3j) ceases to contribute to the 
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Figure 5.1: Evolution of a(rj)H(rj) in a universe with scale factor given by 

Eq. dsm>. 
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power spectrum since the wave lengths of the corresponding GWs exceed 
the size of the horizon. Finally, at 77/1, defined by a (77/1) if (77/1) = a\Hi, all 
GWs have their wave length longer than the Hubble radius and the power 
spectrum vanishes altogether. 

5.2 GWs entropy and its evolution in the dark 
energy era 

There are different expressions in the literature for the entropy density of 
gravitational waves -see e.g., [HHJ EH EH- All of them are based on the 
assumption that the gravitational entropy is associated with the amount of 
GWs inside the horizon. In this section, we review the expressions discussed 
in |Sni- Next, we adopt the proposal of Nesteruk and Ottewill [60 to find 
the evolution to the GWs entropy during the dark energy era. 

5.2.1 Entropy of the GWs 

Branderberger et al. defined the nonequilibrium entropy of cosmological 
perturbations in two ways |55] . One of them is based in the microcanonical 
ensemble jnij while the other is a formula for the entropy that can be as- 
sociated with the stochastic distribution describing the state of the classical 
gravitational field. Both descriptions are in agreement. 

We shall now summarize the second approach concerning GWs. When 
the metric perturbations are sufficiently large, the GWs field can be treated 
as a classical field. The entropy source in this case is the lack of information 
about the exact field configuration. If a classical field ip(t, x) and its canonical 
momentum ir(t, x) are known at all points x at the instant t, the state of 
the system is completely specified and consequently its entropy is zero. But, 
if all we know about the system is the probability of finding it in a state 
(p(t,x), 7r(t,x), i.e., P(<f(t, x), n(t, x)), the entropy can be expressed as 

S = - J ^(V?(x),7r(x)) lnP(^(x),vr(x)) Zty(x)£>7r(x), (5.4) 

where Dip(x.)Dir(x) denotes the functional integral measure for a scalar field. 

When the initial state of the system is Gaussian and the time evolution 
preserves its Gaussian character, the probability distribution P(<p(t, x), n(t, x)) 
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can be expressed in terms of the two point correlations functions {<p(t, x)<^(t, y)), 
{ir(t, x))7r(t, y)) and (f(t, x)7r(t, y)), where (q) stands for the ensemble aver- 
age of q (wich coincides with the space average of q for a spatially homoge- 
neous stochastic process). In this the case, the entropy per unit volume can 
be expressed as 

^ 3 kln((|^| 2 )(|7r fc | 2 )-(|^| 2 |7r fe | 2 )), (5.5) 

where ipk and are the Fourier transformed of <p(t,x) and vr(t,x) respec- 
tively. 

The correlation functions of the GWs field can be expressed in terms of 
the Bogoliubov coefficients (see for more details). Finally, the entropy 
density reads 

d 3 k\nN(k), (5.6) 



where N(k) = | /3Tr2 1 2 is the number of GWs created from the initial vacuum 
with a given wave number k. The above expression is valid in the classical 
limit, i.e., N(k) > 1. 

In (Hill, Nesteruk and Ottewill proposed a definition of the entropy of 
GWs based on the idea that the number of generated GWs describes the 
capacity of the gravitational field to create matter and is associated with the 
gravitational entropy. Their proposal assume the entropy density is propor- 
tional to the GWs number density, i.e., 

s g = An g , (5.7) 

where n g is the number density of gravitational waves, and A is an unknown 
positive-definite constant of proportionality. We shall now make use of this 
definition to evaluate the GWs entropy and to study its evolution in the 
scenario of the previous section. 

5.2.2 GWs entropy in the dark energy era 

We are interested in the evolution of n g during the dark energy era. The 
number density of GWs created from the initial vacuum state is 



dn g (rj) 



|2 PHv)) 
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where the term in square brackets is the density of states and, as we have 
seen, |/?Tr2| 2 is the number of GWs created. We can obtain n g {rj) by inserting 
Eq. ()5.3|) into Eq. ()5.8|) and integrating over the frequency. 

At the beginning of the dark energy era the GWs number density is 
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768vr 5 V a 2 
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, (5.9) 



where 
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- 1 



is the number density at the transition radiation era-dust era. 

i 

, one has 



For rj < 77/0, i.e. < 
n g (V3 < V < Vfo) = n g (r] 2 ) 
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(5.10) 



(5.11) 



The density number given by Eq. (J5.11|) decreases with the scale factor 

because of two effects: (i) the evolution of the volume considered, which 

grows as a 3 , and (ii) the exit of those waves whose wave length becomes 

longer than the Hubble radius, which appears in the term in square brack- 

_i 

ets. As approaches ( — ) 2 , this term tends to zero. From this instant 
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(defined as 77/0) on, the number density reads 
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As time goes on, the term in brackets tends to zero and at the instant rjfi, 
n g vanishes. 

Consequently, the entropy density, proportional to the number density, 
decreases during the dark energy era not just because of the variation of the 
volume considered but also because of the disappearance of the GWs from 
the Hubble volume. 

The GWs entropy inside the event horizon is 

Sg = -^RjjSg , (5.13) 

where Rh = a(t) f°°dt'/a(t') is the radius of the event horizon and t the 
cosmic time. For the horizon to exist Rh must not diverge. For expansions 
of the general form a(rj) oc rf where / < — 1, the horizon exists and can be 
expressed as Rh = —IH^irj). Since it is of the order of magnitude of the 
Hubble radius, H" 1 , we will neglect — I and use both terms interchangeably. 
Making use of Eq. (JE2J) and Eqs. (f5~TT]) - ([5~T2|) , we obtain 
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and 



S g(Vfo < V<Vfi) 
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The GWs entropy is a decreasing function of the scale factor and conse- 
quently of conformal time. The next section explores whether this entropy 
descent can be compensated by an increase of the entropy of the other con- 
tributors, namely, matter and horizon. 
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5.3 The generalized second law 

In this section, we study the implications the GSL of thermodynamics has 
over the GWs entropy during the dark energy era. 

By extending the GSL of black hole spacetimes [021 to cosmological set- 
tings, several authors have considered the interplay between ordinary entropy 
and the entropy associated to cosmic event horizons According the 

GSL of gravitational thermodynamics the entropy of the horizon plus its sur- 
roundings (in our case, the entropy in the volume enclosed by the horizon) 
cannot decrease. Consequently, we must evaluate the total entropy to see its 
evolution during the present dark energy era. 

The cosmological event horizon has an associated entropy that may be 
interpreted analogously to the entropy of the black hole [S5]- An observer 
living in an accelerated FRW universe has a lack of information about the 
regions outside its event horizon (see figure l5~2~j) . The area of the horizon, 
A = Air H~ 2 , represents a measure of this lack of information. As in the case 
of the black hole horizon, we define the entropy of the cosmological horizon 
as Sh = A/A. 

Thus, the entropy of the horizon in our case reads 



and increases with expansion. (Bear in mind that / < — 1, i.e., the dark 
energy behind the acceleration is not of "phantom" type [OS]). 

Here we assume that the entropy of the dark energy field responsible for 
the acceleration does not vary. 

The non-relativistic matter fluid entropy must be also considered. If the 
latter consists in particles of mass m and each of them contributes ks — 1 
to the matter entropy, we get 



for the entropy of the non-relativistic fluid. Here, we made use of the conser- 



vation equation p m {rf) = p m (r) 3 ) ( ) = ^fff ( J with p m the energy 



density of matter. From (|5.17jl . it is apparent that S m decreases with expan- 
sion. 




(5.16) 




(5.17) 
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P 



Figure 5.2: The event horizon of an observer in an accelerated FRW universe 
(a De Sitter universe in this case). The signal emitted by a particle P at the 
instant t < t\ arrives to the observer at the finite instant t*. Therefore, as 
t — > ti, it follows that £* — > oo and the light signal takes longer and longer to 
arrive from the particle to the observer. Signals emitted once the particle is 
beyond the horizon will never reach the observers position. (Adapted from 
Fig. 23.5 of Ref. |53). 
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In virtue of the above equations, the GSL, S' m + S' g + S' H ^ 0, where the 
prime indicates derivation with respect to r], can be written as 
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For / < — 1, both conditions are of the type / [j^-J ^ 0, / being an 

increasing function of a(r]). Therefore, if the condition holds true at the 
beginning of the dark energy era, r\ — 773, it will hold for r) > 773. 

By setting a(r]) = 03 in Eq. (j5.18j) . a restriction over the unknown con- 
stant A of proportionality follows 
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implying that for the GSL to be satisfied the above upper bound must be met. 
In this case the event horizon soon comes to dominate the total entropy and 
steadily augments with expansion. So, even though the entropy of matter 
and GWs within the horizon decrease during the present dark energy era, 
the GSL is preserved provided Eq. (|5.2()jl holds. Note that since the authors 
of Ref. [60J left the constant A unspecified restriction ()5.2()j) turns to be all 
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the more important: it is the only knowledge we have about how large A 
may be. 

Obviously, our conclusions hang on the expression adopted for the entropy 
density of the gravitational waves. Here we have chosen ()5.7j) since, on the one 
hand, it is the simplest one based on particle production in curved spacetimes 
|22j . and on the other hand, s g cannot fail to be an increasing function of n g . 
We believe, that any sensible expression for s g should not run into conflict 
with the GSL, and that the latter may impose restrictions on the parameters 
entering the former. 

As mentioned above, we have left aside models of late acceleration driven 
by dark energy of "phantom type" (i.e., — 1 < I < 0) [HE!- In this case, owing 
to the fact that the dominant energy condition is violated, the event horizon 
decreases with expansion. We shall focus our attention on this issue in the 
next chapter. 

5.4 Conclusions 

The number of gravitational waves inside the cosmological event horizon in 
a four-stage model with scale factor given by (j5.1J) is a decreasing function 
of time in the present dark energy era. The GWs entropy in the horizon 
is associated with the amount of GWs inside it. Consequently, the GWs 
entropy in the horizon decreases with time. 

In order to study the GSL in this scenario, we must consider all the 
possible contributions to the total entropy: the entropy of the horizon itself 
(which is proportional to the horizon area and, consequently, it increases 
with time), the entropy of the dark energy field (which is zero as we assume 
the field is in a specified state) and the entropy of ordinary matter (which 
is proportional to the number of particles of ordinary matter present inside 
the horizon and is a decreasing function of time). If we assume the relation 
between the GWs entropy and the number of GWs is ()5.7|) . the GSL imposes 
a bound over the constant A of proportionality. 
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Chapter 6 



The generalized second law in 
universes dominated by dark 
energy 

In this chapter we leave the main subject of the thesis, the spectrum of 
GWs, to explore the thermodynamics of dark energy taking into account the 
existence of the observers event horizon in accelerated universes. Here we 
shall not consider GWs as we will focus our attention on one single stage 
of expansion. Thus, there is no transitions between different stages and no 
GWs creation. 

As we shall see, phantom and non-phantom dark energies satisfy the gen- 
eralized second law of gravitational thermodynamics whether the equation 
of state parameter is constant or not. Phantom energy transcends the holo- 
graphic principle, non-phantom energy complies with it [68 . 

6.1 Phantom dark energy 

As mentioned before, nowadays there is an ample consensus on the observa- 
tional side that the Universe is undergoing an accelerated expansion likely 
driven by some unknown fluid, the DE, with the distinguishing feature of 
violating the strong energy condition, p + 3p > EH1 E3 EH]. The 
strength of this acceleration is presently a matter of debate mainly because 
it depends on the theoretical model employed when interpreting the data. 
While a model independent analysis suggests it to be below the De Sitter 
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value jZO| it is certainly true that the body of observational data allows for a 
wide parameter space compatible with an acceleration value larger than De 
Sitters jM] • If eventually this proves to be the case, the fluid driving the ex- 
pansion would violate not only the strong energy condition but the dominant 
energy condition, p + p > 0, as well. In general, fluids of such characteristics, 
dubbed phantom fluids, face some theoretical difficulties on their own such 
as quantum instabilities [71]. Nevertheless, phantom models have attracted 
much interest, among other things because in their barest formulation they 
predict the Universe to end from an infinite expansion in a finite time (big 
rip), rather than in a big crunch, preceded by the ripping apart of all bound 
systems, from galaxy clusters down to atomic nuclei jHEj- While this scenario 
might look weird, given our incomplete understanding of the physics below 
p + p = and the scarcity of reliable observational data, it should not be 
discarded right away; on the contrary we believe it warrants some further 
consideration. 

Recently, it has been demonstrated that if the expansion of the Universe 
is dominated by phantom fluid, black holes will decrease their mass and 
eventually disappear altogether [72]. At first sight this means a threat for 
the second law of thermodynamics as these collapsed objects are the most 
entropic entities of our world. This short consideration spurs us to explore 
the thermodynamic consequences of phantom-dominated universes. In doing 
so one must take into account that ever accelerating universes have a future 
event horizon (or cosmological horizon). Since the horizon implies a classi- 
cally unsurmountable barrier to our ability to see what lies beyond it, it is 
natural to attach an entropy to the horizon size (i.e., to its area) for it is a 
measure of our ignorance about what is going on in the other side. However, 
this has been proved in a rigorous manner for the De Sitter horizon only [65 
which, in addition, has a temperature proportional to the Hubble expansion 
rate. Nevertheless, following previous authors [531 EH EH] here we conjecture 
that this also holds true for non-stationary event horizons. 

Phantom expansions of pole-like type 

a(t) oc — (t < t*, < n = constant) , (6.1) 

[t* - t) n 

as proposed by Caldwell arise when the equation of state parameter 
w = p/p = 7 — 1 = constant < — 1. Current cosmological observations hint 
that w may be as low as —1.5 |BTJ] W7\ . In what follows it will be useful to 
bear in mind that n = — 2/[3(l + w)] > 0. 
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From the above equation for the scale factor it is seen that the Hubble 
expansion rate augments 

a n , 

H v = - a = —f (6 - 2) 

whereby the radius of the observer's event horizon 



R H = / ~~ 7T77 = O(t) / — r = f — = i/" 1 (c = 1) 



«(^') Jt a (*') 1 + " 1 + " 



(6.3) 



decreases with time, i.e., Rh < 0, and vanishes altogether at the big rip time 
£*. Consequently the horizon entropy, 

S H = ^ = ttR 2 h , (6.4) 

diminishes with time, Sh < 0. This is only natural since for spatially flat 
FRW phantom-dominated universes one has 

H = -4tt (p + p) > 0, (6.5) 



as these fluids violate the dominant energy condition. Then the question 
arises, "will the generalized second law (GSL) of gravitational thermody- 
namics, S + Sh > 0, be respected?" According to the GSL the entropy of 
matter plus fields inside the horizon plus the entropy of the event horizon 
cannot decrease with time. For the interplay between ordinary matter (and 
radiation) and the cosmological event horizon, see Refs. j^Zl HUB EH HH] • 

6.2 Phantom and non-phantom dark energy 
with constant w 



In this section we study the evolution of the phantom DE with constant w, 
next we compare it with the entropy of the non-phantom DE with constant 
w and, finally, we examine the GSL in both models. 
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In the previous chapter we assumed that the dark energy has no entropy 
associated with it, which would be a reasonable choice if the dark energy is 
given by a scalar field in a pure state. But this is far from being the more 
natural assumption. Henceforward, we shall see both phantom and non- 
phantom fluids as phenomenological representations of a mixture of fields, 
each of which may or may not be in a pure state but the overall (or effective) 
field is certainly in a mixed state and therefore entitled to have an entropy. 
This is the case, for instance, of assisted inflation [71] . 

The entropy of the phantom fluid inside the cosmological event horizon of 
a comoving observer can be related to its energy and pressure in the horizon 
by Gibbs equation 



TdS = dE + pd(±%R 3 H ) , (6.6) 

where owing to the fact that the number of phantom particles inside the 
horizon is not conserved we have set the chemical potential to zero. From 
the relation E = (4ir/3)pR^ I , together with the Friedmann equation p = 
(3/8tt)H 2 , Eq. ()6.3|) . and the equation of state p = wp (with w < —1), we 
get 

TdS= dR H . (6.7) 

1 + n 



Since dRn < the phantom entropy increases with expansion (so long as 
T > 0). 

To proceed further, we must specify the temperature of the phantom 
fluid. The only temperature scale we have at our disposal is the temperature 
of the event horizon, which we assume to be given by its De Sitter expression 

though in our case H ^ 0. Thus, it is natural to suppose that T oc Th 
and then figure out the proportionality constant. As we shall see below, this 
choice is backed by the realization that it is in keeping with the holographic 
principle [ZSj- The simplest choice is to take the proportionality constant 
as unity which means thermal equilibrium with the event horizon, T = Tjj. 
Otherwise energy would spontaneously flow between the horizon and the 
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fluid (or viceversa), something at variance with the FRW geometry. Note 
that, explicitly or implicitly, every FRW model (including inflation models) 
assumes thermal equilibrium. 

After integration of Eq. ([6.7)1 . bearing in mind that S — > as t — > t*, the 
entropy of the phantom fluid can be written as S = —A/ '4. 

Some consequences follow: (z) The phantom entropy is a negative quan- 
tity, something already noted by other authors [TH], and equals to minus 
the entropy of a black hole of radius Rh- (ii) It bears no explicit depen- 
dence on w and, but for the sign, it exactly coincides with the entropy of 
the cosmological event horizon, (m) Since —S oc E 2 it is not an extensive 
quantity. Note that S does not scale with the volume of the horizon but 
with its area. One may argue that the first consequence was to be expected 
since (for T > 0) it readily follows from Euler's equation T s = p + P, s 
being the entropy density. However, it is very doubtful that Euler's equation 
holds for phantom fluids since it is based on the extensive character of the 
entropy of the system under consideration [77] and we have just argued that 
for phantom fluids this is not the case. The fact that S < tells us that 
the entropy of the phantom fluid is not to be understood as a measure of 
the number of microstates associated to the macroscopic state of the fluid 
and that, in general, the current formulation of statistical mechanics might 
not apply to systems that violate the dominant energy condition. Before 
closing this paragraph, it is noteworthy that negative entropies are not so 
new after all; in principle monoatomic ideal gases may exhibit such a feature 
for sufficiently low temperatures [75] . 

Two further consequences are: (iv) The sum S + Sh vanishes at any time, 
therefore the GSL is not violated -the increase of the (negative) phantom 
entropy exactly offsets the entropy decline of the event horizon, (v) | S \ 
saturates the bound imposed by the holographic principle [73]. The latter 
asserts that the total entropy of a system enclosed by a spherical surface 
cannot be larger than a quarter of the area of that surface measured in 
Planck units. (For papers dealing with the holographic principle in relation 
with dark energy, see [7H] and references therein). In this connection, it is 
interesting to see that if the equation of state were such the entropy obeyed 
S* = xS, with k a positive-definite constant, then the GSL would impose 
1 < h. This leads us to conjecture that the entropy of phantom energy is 
not bounded from below but from above, being — Sh its upper limit. 

Nevertheless, in a less idealized cosmology one should consider the pres- 
ence of other forms of energy, in particular of black holes and the decrease in 
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entropy of these objects (Z2j- It is unclear whether in such scenario the GSL 
would still hold its ground. Nonetheless, one may take the view that the GSL 
may impose an upper bound to the entropy stored in the form of black holes. 
At any rate, the calculation would be much more involved and we leave it 
to a future research. Among other things, one should take into account that 
the scale factor would not obey such a simple expansion law as and 
that the black holes would be evaporating via Hawking radiation [32] which 
would also modify the expansion rate and, accordingly, the horizon size. 

A straightforward and parallel study for a non-phantom dark energy- 
dominated universe with constant parameter of state (lying in the range 
— 1 < w < —1/3), shows that in this case too S + Sh = 0. There are two 
main differences, however; on the one hand the entropy of the fluid decreases 
while the area of event horizon augments, and on the other hand S > 0. The 
latter comes from adopting the view that the dark energy must vanish for 
T — > (Planck's statement of the third law of thermodynamics) and realizing 
that this happens for t — > oo. Again, S saturates the bound imposed by the 
holographic principle. 

Another important difference between phantom and non-phantom dark- 
energy expansions is that while the former gets hotter and hotter the latter 
gets colder and colder. This implies the following. As the phantom fluid heats 
up its ability to produce black holes via quantum tunnelling may augment 
and the probability for this to occur becomes high near the Planck limit pTTTj . 
This is a very strongly irreversible process since the variation of entropy 
is enormous going from negative values (phantom fluid) to positive values 
(black holes). Obviously, the black hole production would severely alter the 
dynamics of the last stage of expansion. Two competing effects come into 
play, namely: the tendency to produce black holes via tunnelling, and the 
tendency to destroy these black holes via the Babichev effect [72] (also the 
tendency of the black holes to evaporate may play a non-negligible role). 
Depending on which one of these eventually comes to dominate two rather 
different scenarios arise: (i) If the Babichev effect dominates, the play will 
be rehearsed again and again, probably never reaching the big rip (t = £*) 
but the Universe will always be accelerating, (ii) If the quantum tunnelling 
process takes the upper hand, the Universe will exit the phantom regime and 
come back to the old Einstein-De Sitter expansion during the time taken by 
the black hole to evaporate into radiation. This achieved, the universe will 
start a new radiation dominated era. As a consequence the analysis of above 
will be found inconsistent as a whole because if the Universe is to resume 
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a non-accelerating regime no cosmological event horizon would have ever 
existed. 

Coming back to the phantom energy discussion, one may adopt the point 
of view that the phantom temperature is negative [HO]- But, as mentioned 
above, this would destroy the FRW geometry. On the other hand, negative 
temperatures are linked to condensed matter systems whose energy spectrum 
is bounded from above and may therefore exhibit the phenomenon of popu- 
lation inversion, i.e., their upper energy states can be found more populated 
than their lower energy states [ST] . However, all models of phantom energy 
proposed so far consist in some or other type of scalar field with no upper 
bound on their energy spectrum. In addition, while population inversion is 
a rather transient phenomenon the phantom regime is supposed to last for 
cosmological times. Moreover, bearing in mind that w < — 1 and dRn < 0, 
it follows from Eq. ()6.7|) that if T were negative, then the phantom entropy 
would decrease with expansion and the GSL would be violated. 



6.3 Phantom and non— phantom dark energy 
with variable w 

One may argue that the result S+Sh = critically depends on the particular 
choice of the equation of state parameter w to the point that if it were not 
constant, then the GSL could be violated. Here we shall explore this issue 
taking the phantom model with w dependent on time of Sami and Toporensky 

IBS- 

In this model a scalar field, 4>(t), with negative kinetic energy, minimally 
coupled to gravity sources with energy density and pressure given by 

respectively, is adopted. The equation of state parameter 

p ^ + V(</>) 
w = - = — , 

P ±M - V(<f>) 



is lower than —1 as long as p is positive and depends on time. 
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When the kinetic energy term is subdominant ("slow climb" |22j) the 
evolution equation of the phantom field simplifies to 

mm * ^M. (6.9) 

Assuming also the power law potential V(4>) = Vo<j) a with a a constant, 
restricted to the interval < a < 4 to evade the big rip, the equation of 
state parameter reduces to w ~ — 1 — where x = 2n <ft 2 is a dimensionless 
variable, and the field is an ever increasing function of time, namely, 
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(0 < a < 4), (6.10) 
(a = 4), (6.11) 



where the subscript i stands for the time at which the phantom energy begins 
to dominate the expansion. Thus w increases with expansion up to the 
asymptotic value —1. Likewise the scale function obeys 



a{(f) = a,i exp 
As it should, the Hubble factor 
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augments with expansion while the cosmological event horizon 



R H = x^r(i-^ 



X 



H 



(6.14) 



decreases monotonically to vanish asymptotically -see Figure l6~Tl Here T \^-r L , x) 
is the incomplete Gamma function 



As in the preceding section, the entropy of this phantom fluid can be 
obtained by using Gibbs's equation and T = T# = H/2n. Again, we get a 
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Figure 6.1: Radius of the event horizon vs x for a — 1, 2, 3 and 4 from top 
to bottom. In all the cases Rh — * as t — > oo. For convenience sake we have 
set Xi = 1. 
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negative quantity that increases with time 



here C = v / 3/{2V ) (47r)-2+f a -«/ 4 . 

From the above expressions, it follows that 
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Since the quantity in square brackets is positive-definite for any finite x and 
x > 0, the GSL in the form S + Sh > is satisfied. The equality sign occurs 
just for t — ► oo, i.e., when R H vanishes. 

As Fig. 16.21 illustrates, the holographic principle is again transcended in 
the sense that — oo < S < —Sh- 



As an example of non-phantom dark energy with variable w we con- 
sider the Chaplygin gas [HI]. This fluid is characterized by the equation 
of state p = —A/ p. Integration of the energy conservation equation yields 
p = \J A + (5/a 6 ), whereby w = —A a 6 /(A a 6 + B). Here a is the scale 
factor normalized to its present value and A and B denote positive-definite 
constants. Accordingly, this fluid interpolates between cold matter, for large 
redshift, and a cosmological constant, for small redshift. This has led to 
propose it as a candidate to unify dark matter and dark energy |84~j . 

The radius of the cosmological event horizon reads 
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(6.17) 



where 2-F1 is the hypergeometric function and x = Aa 6 /B is a dimensionless 
variable. As figure IB~3l illustrates. Lim x ^ 00 .Rtf = a/3/(87t) A -1 / 4 which is 
nothing but the value taken by if -1 in that limit. This was to be expected: 
for large scale factor the Chaplygin gas expansion goes over De Sitter's. 



The horizon temperature is 
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Figure 6.2: Evolution of the ratio between the absolute value of the phantom 
entropy and the horizon entropy for the model of Ref. [S2l- The graphs 
correspond successively to a values 1, 2, 3, and 4 from bottom to top. For 
convenience sake we have set x« = 1 in drawing the graphs. All the curves 
asymptote to x — 1. 



78 CHAPTER 6. THE GSL IN UNIVERSES DOMINATED BY DE 



1.002 




0.988- 

10 20 30 40 50 60 70 80 

X 



Figure 6.3: Evolution of the adimensionalized radius of the event horizon in 
the Chaplygin gas model. It tends to unity for x — > oo. 
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while the entropy of the fluid and the entropy of the horizon evolve as 
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respectively. The two latter equation imply that 5 + S# > for x > 2.509 
or equivalently, for pfu\ < LISA 1 / 2 . Thus there is an early period in which 
the GSL seems to fail. However, this may be seen not as a failure of the 
GSL but as an upper bound on the initial value of the energy density of the 
Chaplygin gas. This is somewhat similar to the bound imposed on the initial 
entropy of relic gravitational waves in a universe dominated by non-phantom 
dark energy [57] . 

Notice that, S is a positive, ever-decreasing function of the scale factor 
that complies with the holographic bound S < Sh- 



6.4 Quasi- duality between phantom and non- 
phantom thermodynamics 



The duality transformation 

H 

P + P 

leaves the Einstein's equations 

H 2 



H=-H 

P + P= ~(P + P) 



(6.21; 



p + 3H(p'+p) = 0, 
H = -4vr(p + p), 



(6.22) 
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of spatial FRW universes invariant and the scale factor transforms as a — > 
a = 1/a [HH]. Thus, a phantom DE dominated universe with scale factor 
a = l/(— t) n (for simplicity we assume a phantom scale factor as Eq. ()6.1|) 
with t* = and t < that tends to - ) becomes a contracting universe 
whose scale factor obeys a = (—t) n . This universe begins contracting from 
an infinite expansion at the infinite past and collapses to a vanishing scale 
factor for t = 0~ (left panel of Fig. I6.4jl . 

Reversing the direction of time (i.e., the future of t, t + dt, transforms into 
t—dt), the contracting universe is mapped into an expanding one whose scale 
factor grows with no bound as t tends to — oo. Notice that the first trans- 
formation takes the universe from expanding to contracting and the second 
transformation makes the universe expand again, i.e., the final cosmology has 
H > 0. It should be remarked that the latter transformation also preserve 
Eqs. flPH) . 

Under this two successive operations, the final and the initial equation of 
states parameters are related by w — > w = —(w + 2). Consequently, phantom 
DE dominated universes with n > 1, (i.e., —5/3 < w < —1) are mapped into 
non-phantom DE universes (see Fig. 16. 4j) . Therefore, both of them, the 
original and the transformed universes, have an event horizon and entropies 
whose relation we derive below. 

The radius of the event horizon of the original universe is Rh = (—t) /(l + 
n) meanwhile the radius of the event horizon of the transformed universe is 
Rh = (—t)/(n — 1) = Rh[(1 + n )/(l ~ n )]- The transformations preserve 
the horizon temperature as the Hubble factor does not change. From the 
definition of the horizon, it is readily seen that 



In virtue of the above expression together with the Gibbs equation one ob- 
tains 



dRtf 



HR H - 1. 



(6.23) 



dt 



T— = -47r(p + p)R 2 H , 
whereby, the entropy transforms as jHHl 

S -»• S = -[(1 + 3iu)/(5 + 3w)} 2 S. 



(6.24) 



(6.25) 



Therefore, the entropy of the final non-phantom dark energy dominated uni- 
verse is proportional, but of the opposite the sign, to the entropy of the 
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t t+dt t t+dt 



. > < . 

Figure 6.4: The left panel depicts the scale factor of the original phantom 
dominated universe (solid line). The dashed line shows the evolution of the 
scale factor dual of this phantom DE model, a contracting universe. In both 
models (the phantom and its dual) the time grows from left to right as tends 
to t = 0~ and the instant t + dt lies at the future of the instant t. In the right 
panel, the time direction is reversed (t + dt lies at the past of t). Under this 
operation, the contracting model of the left panel becomes the non-phantom 
DE dominated universe of the right panel if —5/3 < w < — 1. 

original phantom universe. We might say that the duality transformation 
"quasi" preserves the thermodynamics of dark energy. This result is in keep- 
ing with the findings of Section 16.21 

6.5 Conclusions 

We explored some thermodynamical consequences for the dark energy when 
the entropy of the cosmological event horizon is not ignored. Irrespective of 
if w is constant or not, we have found that phantom fluids possess negative 
entropy, transcend the holographic bound and their temperature and entropy 



82 



CHAPTER 6. THE GSL IN UNIVERSES DOMINATED BY DE 



increase as the Universe expands. By contrast, non-phantom fluids have 
positive entropy, satisfy the holographic bound and their temperature and 
entropy decrease with expansion. In all the cases, the GSL is fulfilled. 

It goes without saying that negative entropies are hard to assimilate; espe- 
cially because the Einstein-Boltzmann interpretation of entropy as a measure 
of the probability breaks down under such circumstances. Systems of neg- 
ative entropy appear to lie outside the province of statistical mechanics as 
is currently formulated [ED EH] ■ Nevertheless, the laws of thermodynamics 
alone do not entail that S must be a positive quantity. The latter will be 
found to be positive or negative only after combining these laws with the 
equation (or equations) of state of the system under consideration. What 
we believe to be at the core of thermodynamics is the law that forbids the 
entropy of isolated systems to decrease, this law contains the GSL for gravi- 
tating systems with a horizon. 

As pointed out by Nojiri and Odintsov, it may well happen that in the last 
stages of phantom dominated universes the scalar curvature grows enough 
for quantum effects to play a non-negligible role before the big rip. If so, 
the latter may be evaded or at least softened [HZj. While our study does not 
incorporate such effects they should not essentially alter our conclusions so 
long as the cosmic horizon persists. 

It is noteworthy that the duality transformation p + p — > — (p + p) and 
H — > —H along with reversing the direction of time (i.e., dt — > — dt) leaves 
Einstein's equations invariant and maps phantom cosmologies, with both 
H and H positive, into non-phantom cosmologies with H > but H < 
-see, e.g., jHS]. However, this duality does not necessarily extend to the 
thermodynamics of the respective universes since future event horizons exist 
only when a > 0. Nevertheless, in the particular case of pole-like expansions, 
Eq. (|6.ip . the transformation preserves the temperature while the entropy 
transforms as S — » — [(1 + 3w)/(5 + 3w)] 2 S with the constraint —5/3 < w < 
-1. 



Summary 



Cosmology has for a long time been a rather speculative science. Hubble's 
discovery that the Universe is expanding, and -more recently- the realization 
that at present this expansion is accelerated, the measured abundance of 
light elements, the mass distribution of galaxies and clusters thereof, and the 
discovery and posterior measurements of the anisotropies of the CMB have 
changed this picture. Hopefully, measurements of GWs will soon be added 
to this short list. At any rate, now we can speak confidently of physical 
cosmology as a fully-fledged branch of Science. 

The relic GWs constitute a privileged window to determine the evolution 
of the Universe. Little is known from the early evolution of the Universe 
and the predictions for their spectrum depend on the model considered. Ac- 
cording to these predictions, a spectrum of relic GWs is generated making 
feasible its detection with the technology currently being developed. 

In this thesis, using the adiabatic vacuum approximation, we have re- 
viewed how the expansion of the Universe amplifies the quantum vacuum 
fluctuations, and how the relic GWs spectrum is related with the scale fac- 
tor. 

We have later evaluated the spectrum in a four-stage model (which consist 
on a De Sitter stage, a stage dominated by a mixture of MBHs and radia- 
tion, a radiation dominated stage and finally a non-relativistic matter (dust) 
dominated stage). We have demonstrated that the spectrum in this scenario 
is much lower than the predicted by three-stage model (De Sitter-radiation 
era-dust era). We have also shown how the bound over the GWs spectrum 
from the measured CMB anisotropies places severe constraints over the free 
parameters of the four-stage model. 

We have also considered a scenario featuring an accelerated expanding 
era dominated by dark energy, right after the dust era of the three-stage 
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model. We have found that the current power spectrum of this four-stage 
scenario exactly coincides with that of the three-stage, but it evolves in a 
different fashion. We have considered as well the possibility that the dark 
energy decays in non-relativistic matter leading to a second dust era in the far 
future and obtained the power spectrum of the GWs as well as the evolution 
of the density parameter. 

We have applied the generalized second law of thermodynamics to the 
four-stage model of above. Assuming the GWs entropy proportional to the 
number of GWs, we have found the GSL is fulfilled provided a certain pro- 
portionality constant does not exceed a given upper bound. 

Finally, we have extended the GSL study to a single stage universe model 
dominated by dark energy (either phantom or not), and found that the GSL 
is satisfied and that the entropy of the phantom fluid is negative. Likewise, we 
have found a transformation between phantom and non-phantom scenarios 
preserving the Einstein field equations that entails a "quasi" duality between 
the thermodynamics of both scenarios. 



Appendix A 

Sudden transition 
approximation 



Let us assume a spatially flat FRW universe which gradually evolves from 
an initial stage (whose expansion is dominated by a barotropic fluid of index 
7(t—i)) to a next stage dominated by another perfect fluid of barotropic index 
7 r . If the transition between both stages begins at the instant rf = rji — Ar]i 
and ends at 77 + = rji + A^, the scale factor as a function of time in this 
scenario reads 



a(rj) = { a tr (v), (v < V < V + ,) , 

{^T^-^ + ^y, (V + <V), 

(A.l) 

where the subindex i — 1 denotes the instant when the expansions begins and 
the subindex + denotes the quantities evaluated at r/ + . The function a tr {ji) 
describes the smooth transition from one stage to the other. Consequently, 
in this scenario a{rf) and a!{rj) are continuous functions of time at the points 
77 _ and rj + . 

The solution to equation (|2.10j) is 

( v ^/2)e^(^ 1 )A;-V2 x V2 i) ^)_ i) ^ (xi(r _ i))? <v < r ), 

%._!) (v) = \ fr r (r)), (77- < n < r] + , ) , 

aiVir (v) + W r (v) , (V + < V) > 

(A.2) 
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but also 



<^ (r _i)(^) - ft^.^W. (Vi-i <V<V~), 

ViM = { "Ptriv), (V <V < ) , (A.3) 

( y ^/2)e^rk- 1 / 2 x 1 l l 2 H^ 1 _(x lr ), (77+ < V ) , 

where the solutions in the transition region fi tr and JZ tr are some unknown 
functions. 

By imposing the continuity of fii (r _ 1) and fi[ r and its first derivatives at 
rf and r/ + , we get 



an 



= i mLiM^Ut) - vi r -i(v~)~P*tr'(v~) 

Pi = i \a tr (r] + )fi f lr {r] + ) - fi' tr (ri + )fi lr {r] + )] 



(A.4) 



(A.5) 



For modes whose period is much larger than the duration of the transition 
era (i.e., those modes allowed by the adiabatic vacuum bound u{rf) < 0^(77)), 
the transition can be assumed instantaneous. This is the so-called "sudden 
transition approximation". In this limit, we can safely let 77" — > r/ + . Re- 
labelling the instant 77 + as rji, we have that Utrijli) = ^hr-v, (Vi)i "PtriVi) = 
Hi T {r]i), and a similar relation for the derivatives evaluated at r/j. Equations 
(j2.41|) and f)2.42j) readily follow from (|A.4|) and ()A.5|) in this approximation. 
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